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Abstract. The usual treatment of a (first order) classical field theory such 
as electromagnetism has a little drawback: It has a primary constraint sub- 
manifold that arise from the fact that the dynamics is governed by the anti- 
symmetric part of the jet variables. So it is natural to ask if there exists a 
formulation of this kind of field theories which avoids this problem, retaining 
the versatility of the known approach. The following paper deals with a family 
of variational problems, namely, the so called non standard variational prob- 
lems, which intends to capture the data necessary to set up such a formulation 
for field theories; moreover, we will formulate a multisymplectic structure for 
the family of non standard variational problems, and we will relate this with 
the (pre)symplectic structure arising on the space of sections of the bundle of 
fields. In this setting the Dirac theory of constraints will be studied, obtaining 
among other things a novel characterization of the constraint manifold which 
arises in this theory, as generators of an exterior differential system associated 
to the equations of motion and the chosen slicing. Several examples of applica- 
tion of this formalism are discussed: Two of them motivated from the physical 
point of view, that is, electromagnetism and Poisson sigma models, and two 
examples of mathematical application. In the case of electromagnetism, it is 
shown that this formulation avoids the problems arising in the usual approach. 

1. Introduction. In this work, the term field theory will refer to a particular kind 
of variational problems on the sections of some bundle on space-time (for defini- 
tions, see section 2). In the usual approach, the Dirac constraints in a field the- 
ory are found using the Gotay-Nester-Hinds algorithm on an infinite dimensional 
presymplectic manifold associated to the underlying variational problem. One of 
the purposes of this article is to show another way to build up these constraints, 
namely, by using a geometrical representation of the corresponding Euler-Lagrange 
equations into the exterior algebra of a bundle. The tools used for this task were 
taken from the theory of exterior differential systems (EDS from now on, for refer- 
ences see appendix B). In this way the Cauchy-Kowalevsky existence theorem can 
be introduced into the realm of field theory; on the contrary, the usual setting has 
to do with functional analytic conditions, hiding these existence conditions into the 
(infinite dimensional) manifold structure. 
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In more precise terms, the variational problems we are concerned with are initially 
characterized by a double fibration 

A -> Ai -> M n 

together with an exterior differential system I C f2* (A) and an n-form A. These 
kind of problems will be called non standard variational problems throughout the 
work. The idea for this formulation can be traced back to [18], for the case in 
which the base manifold has dimension n ,= 1: another work on the subject, in the 
same vein, can be found in [21]. It is worth remarking that both references provide 
examples from geometry (not only mechanics) where this scheme can be applied. 
The basic idea is that, for every section a of the bundle Ai, we can build a section 
prer for the bundle A by using these data (see details below); the non standard 
problem consists of finding those sections a of the bundle Ai which are extremals 
of the functional 

S x [a] := f (pr*)* A. 

This setting includes what can be called standard variational problems, like the 
usual Lagragian mechanics and the lagrangian viewpoint of first order field theories 
as in [15]: For example, in the latter case the bundle Ai corresponds to the bundle 
of fields, A is the 1-jet space of sections for this bundle, and the prolongation of a 
section a is given by its 1-jet j x o '. 
Our aim is twofold: 

1. In first place, to set up an infinite dimensional presymplectic manifold with 
a hamiltonian for every non standard variational problem. This is done by 
appealing to a bivariant Lepagean equivalent variational problem associated 
to the non standard problem; it provides us with a multisymplectic-like struc- 
ture, and introducing an slicing of the bundle where this structure lives, the 
presymplectic manifold can be defined. 

2. Secondly, we want to prove that the Dirac constraints obtained from the suc- 
cessful application of Gotay-Nester algorithm to the data found in the previous 
item can be calculated as the generators of some EDS closely related to the 
(Lepagean equivalent) variational problem. We must remark here that these 
constraints generalize in some sense the Dirac constraints to the non standard 
setting: The usual Dirac constraints can be found through this procedure by 
using the canonical Lepage equivalent (in the sense of definition 2.4 below) 
for a first order field theory on the 1-jet bundle associated to the bundle of 
fields, as shown in example 11. 

The present article is structured as follows: In section 2 the non standard varia- 
tional problem is defined, and some relevant examples are presented. In particular, 
it is shown that several dynamical problems of mathematical physics can be formu- 
lated as non standard problems; as a particular example, we discuss in this setting 
the electromagnetism. The Euler-Lagrange equations of the non standard problem 
are found here, and the ideas of Lepage equivalent and canonical Lepage equivalent 
system associated to a non standard variational problem are introduced. This gives 
us a kind of multisymplectic structure for every non standard variational problem. 
It is important to note that no use of Legendre transform is made in the construc- 
tion of this multisymplectic space; the main disadvantage of this approach is that 
the equivalence between the original equations of motion and the equations of mo- 
tion in the multisymplectic space must be done in each case separately (this is an 
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issue related to the so called bivariance of the chosen Lepage equivalent problem, 
see below). By adding a compatible slicing of the space-time, a Hamiltonian version 
of the equations governing the extremals is developed in section 3: it is described a 
procedure that associates a presymplectic manifold (infinite dimensional if n > 1) 
and a function on it to every (Lepagean equivalent of a) non standard problem. 
Moreover it is proved here that the solutions of the dynamical system determined 
by these data are extremals of the underlying variational problem. Section 4 con- 
tains the main result of the article, which deals with the description of the constraint 
submanifold arising from the Gotay-Nester algorithm [ LI ] in terms of an exterior 
differential system associated to the data of the theory. The section 5 contains 
four examples where these considerations are applied. The first two constitues the 
examples with physical content, namely, electromagnetism and the Poisson sigma 
model, and the last two explores other aspects of the method: the calculation of 
integrability conditions for a system of PDEs via Gotay-Nester algorithm, and a toy 
model where this algorithm fails to reach a succesful termination. The relevance of 
the last example is that its features could be related to the singular behaviour of 
some field theories (cf. remark 5.7). We show also how the EDS perspective help 
us in the interpretation of the results. A note of caution about these examples is in 
order here: It will be assumed that there exists a well-defined differential structure 
on the sets of sections which we work with, so several of the manipulations in this 
context are of formal nature. 

It is worth noting that the idea giving rise to the definition of the EDS associated 
to the constraints is an extension of the one used in [19] for the treatment of clas- 
sical mechanics from the EDS viewpoint. Similar results, but from the viewpoint 
of Janet-Riquier theory, instead of Cartan's EDSs, and oriented to the construction 
of an algorithm, can be found in [28, 29, 12, 11]. The introductory material about 
variational calculus and exterior differential systems was deferred to Appendices A 
and B. 

Notation. Given a bundle F — > B , T (F) denotes the set of differ entiable sec- 
tions of this bundle. The subbundle VF C TF denotes the set of vertical vectors. 
The symbol f2* (X) indicates the set of forms of any order on the manifold X , and 
correspondingly f2 fe (X) is the set of k-forms on X . Although in some cases this 
fact was proved (see proposition 6 below), it is assumed that some vector fields de- 
fined on certain subsets of the bundles which we work with admits an extension to 
a neighborhood of these subsets. We do not discuss the convergence of the inte- 
grals used elsewhere; it is supposed that the hypothesis ensuring its convergence are 
satisfied. The notation Tx (resp. fx) used in [1] in order to denote the canoni- 
cal projection of the tangent bundle (resp. cotangent bundle) of X is adopted. If 
i : Y °-> X is the canonical infection of some submanifold Y of X , the symbol 
TyX will denotes the pullback bundle i* (TX). Given a set of forms S C fi* (X) 
we denote by (S) alg C f2* (X) the algebraic ideal generated by these forms, and by 
{S) di g C f2* (X) the minimal algebraic ideal containing S and is closed with respect 
to the exterior derivative operator. If S C V is a set of vectors in the vector space 
V , we indicate by (S) the subspace spanned by the vectors in S. 

2. Formal structure for variational problems. The essential data in both the 
formulation of the usual field theories [15, 16, 2] and classical mechanics [1, 19] are 
the following: A double fibration A — ^-s> Ai — ^ M (with composition 7r := tti op), 
an exterior differential system Ic!]' (A) and a semibasic n-form A G VL n (A), where 
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n = dim M. In terms of these objects, the prolongation of a section a £ T (Ai) can 
be defined as a section prcr g T (A) such that 

1. Its graph is an integral submanifold for X, that is, prcr* (I) = 0, and 

2. The section prcr covers cr, that is, the following diagram is commutative 



A 




M -5— Ai 

Hence we introduce the action as a real- valued function on the set of sections of Ai 

S[a}:= [ (pra)* (A). 

Definition 2.1 (Non standard variational problem). In the previous setting, the 
non standard variational problem consists in finding sections cr 6 T (Ai) which are 
extremals of the action S. 

The usual point of view in geometric mechanics is to take the contact structure 
on jet manifolds as the basic prolongation structure for the underlying variational 
problem; in this view it is a piece in the formalism. The non standard approach 
consists in the weakening of this assumption, regarding the prolongation structure 
as an additional data of the problem, related to some aspects of its underlying 
geometry. The following examples may help us to clarify this assertion. Moreover, 
they are in order to illustrate the ubiquity of this notion; subsequent sections in the 
article will intend to show its usefulness. 

Example 1 (Classical Mechanics). Taking M := J C 1 a real interval, Q a man- 
ifold, Ai := Q x I and A := TQ x /, with tt\ (q,t) = t and p(q, q;t) = (q,t), 
choosing as the prolongation EDS the differential ideal generated by the forms 
6 l := dq l — q l dt, and with A := Ldt, we obtain the usual variational problem of the 
classical mechanics. In this case if we have that a : 1 1-> (q (t) ,t), then 

pra : t M. (q (t) , ^ (f) ; tj 

is the formula defining the prolongation of sections. 

Example 2 (First order field theories). The following example show how first order 
field theories fits in our scheme, that is, how one might formulate a first order field 
theory as a non standard problem. So let us suppose that the fields are sections of 
certain bundle F — M on space-time. The double fibration is in this case 

A := J 1 (tt) ^4 A := J° (tt) = F -> M, 

where J k (tt) denotes the fc-jet manifold of F. By taking on J 1 (tt) the contact EDS 
as prolongation structure, it can be shown that prcr = If L £ C°° (J 1 (tt)), 

the n-form A := Ldx 1 A • • • A dx n satisfies 

(pra)* A = (L o fa) dx 1 A • • • A dx n ; 

the variational problem consisting in finding the extremals to 

cr i y I (Lo j x a) dx 1 A • ■ • A dx" 

JM 
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is the usual Hamilton's principle for the first order field theory with lagrangian 
density L. 

Example 3 (Electromagnetism). Let us consider now the first example on non 
standard variational problem, that is, the electromagnetic field on a spacetime M; 
previous work on the subject can be found in [9, 20]. Let us take 

2 

A: = f\(T*M)(BT*M 
Ai : = T*M 

7Ti := f M and for p the projection p 2 : f\ 2 (T*M) © T*M -> T*M onto the second 
summand. An arbitrary element of A is given by a pair (F, A) where F <E A CCM) 
and A £ T^M for some m £ M. For every k € N there exists a canonical /c-form 
on /\ k (T*M) defined by 

k 

@k\ a --=aof M *, Vae/\(T*M); 
as prolongation EDS X we take the differential ideal generated by the form 

T := d6i - e 2 . 

So a section mi4(F (m) , A (m)) of A is an integral section for X iff F = dA. 

Example 4 (EDS as a non standard problem). This is the first example of a non 
standard problem of mathematical nature: Every EDS I on a bundle F — > M can 
be considered as a non standard problem, by using the double fibration 

F M M, 

with prolongation structure X and lagrangian n-form A = 0. Therefore the interme- 
diate bundle has just one section, namely oo = idM, and the possible prolongations 
of it are the integral sections of X. Later on we will use this remark in order to 
employ the Gotay-Nester algorithm in searching the integrability conditions of a 
system of PDEs (see example in section 5.5). 

Example 5 (Tetrad gravity). Let M be a pseudoriemannian 4-manifold. Let us 
introduce a formulation of the variational problem for tetrad gravity on M in these 
terms. The fields in this case are the components e M e 57 1 (M) ,fj, = 0, • ■ ■ ,3 of 
a tetrad on the spacetime M. The natural "velocities" for such an object are not 
its prolongation into J 1 (tm)\ instead, it is more geometric to take as velocities 
for a tedrad the associated connection forms £ fi 1 (M). Thus it is possible to 
formulate tetrad gravity as a non standard problem on the double fibration 

Fa (A/) © (T*M) m -> F a (M) -> M 

where pM '■ Fq (M) — > M stands for the coframe bundle on M, with prolongation 
structure generated by the forms 

9 Q := dE a — Q,™ A E v , u = 0, - ■ ■ ,3. 

Here with capital letters we denote the corresponding canonical forms. The la- 
grangian is the scalar curvature, so 

A := * (E a A E' 3 ) A (dn aP + Vl u a A Sl vp ) . 

The prolongation means that it is a torsionless connection. 
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The purpose of the forthcoming sections is to show that in the non standard 
setting it is possible to build a reasonable multisymplectic space and even to deal 
with its Dirac constraints, althought no use of Legendre transformation will be made 
in the process. But first let us show how the Euler-Lagrange equations for a non 
standard variational problem can be built. 

2.1. Euler-Lagrange equations for non standard variational problems. In 

order to perform variations on the class of functionals relevant to our study, we 
need to adapt the proposition 7 to the case in which the section involved in the 
integrand is the prolongation of another one by using some prolongation structure. 
To this end, let us introduce the following definition. 

Definition 2.2. An (infinitesimal) allowed variation for a non standard varia- 
tional problem with prolongation structure given by the EDS I is an element 
V E r(cr* (FA)), a e T (1) admitting an extension V G T(VA) C X(A) which 
is an infinitesimal symmetry for I, that is 

a* (Cyl) = 0. 

It is assumed, unless explicitly stated, that whenever a manifold has a boundary, 
the allowed variations annihilates on this boundary. Then we have the following 
proposition [13, 18]. 

Proposition 1. The Euler-Lagrange equations associated to the functional S [a] = 
J M (prer) A defined on T (Ai) are given by 

(pra)* (VjdX) = mod exact n-forms on A, 
where V G T ((pre)* (VA)j is an arbitrary allowed variation. 

Example 6 (Electromagnetism - Cont.). We return to the electromagnetic field in 
order to find the Euler-Lagrange equations for the lagrangian A :— 82 A *82; here 

*&2\a ■= (*£*) ° Tm* 

by using the Hodge star associated to the spacetime metric. Now, as in every 
vector bundle, the vertical space can be identified at each point with the fibre 
through that point, so any vertical vector field on A can be denoted as an application 
(F,A) G A h- > (6F,8A) G ir^ 1 (it (F, A)) ~ V {FtA) A. This identification has the 
following interesting property: For the Lie derivative the following holds 

k 

{£ Sa O k )\ a = foof M „ Vae/\(T*M). 

So if V := (6F, 5 A) is an allowed variation with extension V := \5F, 5Aj to a 
vertical vector field on A, this extension must be such that 

£{ S fm) (de 1 -e 2 ) = / x(de 1 -e 2 ) 

for some function fj, 6 C°° (A). By restricting to the integral submanifold Im (prer) , a € 
r (X) we obtain the following differential condition 

d5A -5F = 0. 

Then the Euler-Lagrange equations are expressed as 

2d<5AA*F = 8 A A d (*F) = <(=> d(*F) = 

to which the prolongation condition must be added 

F = dA. 
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2.2. (Multi)hamiltonian formalism through Lepagean equivalent prob- 
lems. We want to construct a hamiltonian version for the non standard variational 
problem. The usual approach [15] seems useless here, because of the following facts: 

• The covariant multimomentum space is a bundle in some sense dual to the 
velocity space, which is a jet space. 

• The dynamics in the multimomentum space is defined through the Legendre 
transform, and it is not easy to generalize to a non standard problem this 
notion. 

The trick to circumvect the difficulties is to mimic the passage from Hamilton's 
principle to Hamilton-Pontryaguin principle. This is done by including the genera- 
tors of the prolongation structure in the lagrangian density by means of a kind of 
Lagrange multipliers. This procedure will be formalized below, where the hamil- 
tonian version is defined by associating a first order variational problem to the non 
standard variational problem, whose extremals are in one to one correspondence 
with its extremals. This is called canonical bivariant Lepage equivalent problem. 

2.2.1. Lepagean equivalent problems. Here we will follow closely the exposition of 
the subject in the article [13]. Before going into details, let us introduce a bit of 
terminology: If A — ^> M is a bundle, A 6 O™ (A) (n = dim M) and X is an EDS 
on A, the symbol ^A — M,X,a\ indicates the variational problem consisting in 
extremize the action 

S[a}= [ a* (A) 

J M 

with a € r (A) restricted to the set of integral sections of X. Furthermore, £ (A) 
will denote the set of extremals for (A M,X, A). 

The idea is to eliminate in some way the constraints imposed by the elements 
of X: intuitively, it is expected that the number of unknown increase when this is 
done. The following concept captures these ingredients formally. 

Definition 2.3 (Lepage equivalent variational problem). A Lepagean equivalent of 
a variational problem ^A — M, X, \ j is another variational problem 

(A^M,{0},A) 

together with a surjective submersion v : A — > A such that 

• p = 7r o z;, and 

• if 7 € F (jCj is such that v o 7 is an integral section of X, then 

7*A = (v o 7)* A. 

There exists a canonical way to build up a Lepage equivalent problem associated 
to a given variational problem ( A — — > A/, X, \ j , the so called canonical Lepage 
equivalent problem. Let X be differentially generated by the sections of a graded 
subbundle I C f\' (T*A) (this is a "constant rank" hypothesis, ensuring the existence 
of a bundle in the construction, see below). Define I alg as the algebraic ideal in 
f2* (A) generated by F (/), and 

(Z alg )' :=Z alg nO' (A). 
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For X e n n (A), define the affine subbundle W x C A" ( r * A ) whose nber above 
p G A is 



W x \ := 



{X\ p +(3\ p :l3e(l^) n }. 



Definition 2.4 (Canonical Lepage equivalent problem). In the previous setting, 
it is the triple (\V X -A M, {0} , 0\ , where v is the canonical projection : 

A™ (T*A) —> A restricted to W x , p := tt o v and is the pullback of the canonical 
n-form 

9»l„:=«°(^), 

to W x . The form 9 will be called Cartan form of the variational problem. 

It is worth remarking here that the terminology adopted here correspond to [13], 
which is slightly different from the terminology of the classical theory, as exposed 
in e.g. [24]; it is fully explained in the former work how can be related the notions 
developed in each case. 

Returning to our main concern, it can be proved that the canonical Lepage equiv- 
alent is a Lepagean equivalent problem of ^A — M,X, \ j . Now, the extremals of 
some variational problem has, in general, nothing to do with the extremals of its 
Lepagean equivalent problem, so it is necessary to introduce the following definition. 

Definition 2.5 (Covariant and contravariant Lepage equivalent problems). We 
say that a Lepagean equivalent problem (w x — ^ M, {0} , Bj for the variational 

problem ^A — M, I, Xj is covariant if v o 7 G £ (A) for all 7 £ £ ; on the 
contrary, it is called contravariant if every a £ £ (A) is the projection of some 
extremal in £ ^Qj through v. A Lepagean equivalent problem is bivariant if and 
only if it is both covariant and contravariant. 

There exists a fundamental relation between the extremals of a variational prob- 
lem and the extremals of its canonical Lepage equivalent. 

Theorem 2.6. The canonical Lepage equivalent is covariant. 

For a proof, see [13]. The contravariant nature of a Lepage equivalent problem 
is more subtle to deal with. Next we will describe another construction for the 
Lepage equivalent problem, by reducing the number of additional variables that we 
put in order to get the prolongation EDS into the variational equations; although 
in this way there exist more chances to find a non contravariant Lepage equivalent 
problem, we will see that in many cases it leads to well-behaved Lepage equivalent 
problems. 

2.2.2. Another canonical Lepage equivalent problem. The canonical Lepage problem 
has to deal with general EDSs; the examples we will try to manage has some nice 
features that allow us to simplify this scheme. Namely, define the subbundle 



Zi (A) := ja e f\ (T*A) : wijtfcja 



Vui,«2 (A — s- M) -vertical vectors 



and suppose further that the generators of X are in Z\ (A); as before, we have a 
subbundle (X al «)" C A™ (T*A). Let K n C I n be the subbundle whose fibers are 
the subspaces K n \^ C I n \ 1 of forms of (algebraic) degree 1 in the generators; set 
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K n := K n n Zi (A). Depending on I, it can leads to an affine subbundle W x in 
A* (T*A) with fibers 



W 



XL + tl> : 4' 6 K n 



As before, the map v : W x — > A is the restriction of the canonical projection 
: /\* (T*A) — > A; moreover, it can be shown that it is a covariant Lepagean 
equivalent problem also. The reasons for this choice will become apparent later: 
The variational problems associated to forms in Z± (A) will have a well defined 
description in terms of dynamics on a (infinite dimensional) presymplectic manifold 
(cf. prop. 3.4). 

A local version of the alternate canonical Lepage equivalent. Let us now define the 
pullback bundle tt* (A* (T*M)) on A; then we have the commutative diagram 



f\ (T*M) 



f\ (T*M) 



A 



M 



Let Z (A) C A* (T*A) be the subbundle of semibasic forms on A, and Z\ (A) := 
Z Q (A) n A fe (T*A) for all k e N. The following lemma can be found in [27] (cf. 
lemma 3.3.5 there). 

Lemma 2.7. There exists a natural identification between the bundles n* (A* (T*M)) 
and Zq (A) . 

This identification provides us the horizontal map /, which has the following nice 
property. 

Proposition 2. The submersive bundle map f : Zq (A) — > A* (T*M) is graded (of 
degree ), and verifies that 

(/l4(A))*ef =e£|z*(A), vfceN (i) 

where 9^ G n k [f\ k (T*A)j and 9f e fi fe (A fc (T*M)^ are tfte canonical k-forms 
on the corresponding spaces. 

Further simplifications can be performed whenever the EDS admits a set of global 
generators 

J l 7 = K ^"1| 7 >■■■ , ajfcj 7 , ai| 7 r- - , afe 2 | 7 >" ' » a il 7 >••' » a fe P 7 ^ ; V l e A > 
where (cej, • • • , oP k . ) = =: Ifl O-' (A). It is supposed here that a{ e Z\ (A) 

1 1 (WW 

1, • • • , ftj, although neither of them is in Zq (A). Let us 



for all j = 1, 



3 / diff 



define A := A Xj/ ©f =1 (A™' T*M) 1 where 1 to; := n — I; the previous lemma 



1 Without loss of generality, it will be assumed that p < dimAf = n, because we are dealing 
with integral sections, a subset of integral manifolds of dimension n for X. 
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leads us to the map F\ : W x — > A given by 

1=1 i=l 

* (7, / (/tj >•••,/ >•••>/ (/t J >•••>/ (/&,)) • 

Let us introduce the definitions 

mi 

Z mi := Z (A) n /\ (T*A) , 2*; := Z m x • • • x Z m; 

ki times 

for I = 1 , ■ • • , p. The map F\ will be well-defined provided that between the gener- 
ators there are no "zyzygies of degree n", namely, that the unique solution for the 
system 

EE«!|/ T ;=o 

1=1 i=l 

in Z^ ± x • • • x Zm p is 7* = 0. Let us suppose that it is indeed the case; on A we 
will define the n-form 

p / fc; \ 

In this setting the proposition 2 implies the following result, that allow us to use 
A as canonical Lepage equivalent whenever our prolongation EDS I admits global 
generators. 

Lemma 2.8. Let us suppose that I admits a set of global generators with the proper- 
ties detailed above; then the map F\ : W x — > A is a diffeomorphism, and F^X = 0. 

These considerations can be applied to build a local version for the canonical 
Lepage equivalent variational problem, because it is always possible to find an open 
set where the considered EDS is described by a set of forms as above, so this pro- 
cedure provides a local description (i.e. on the given open set) for the modified 
canonical Lepage equivalent. 

As we said before, the main difference between our constructions and the construc- 
tion of the canonical Lepage equivalent problem is that in our case the affine bundle 
has lower rank. This means that there are less variables to take care of; the main 
drawback is that such a system has lower chances to be contravariant. Nevertheless, 
it will be shown below that it works very well in many important circunstances. 

2.2.3. Canonical Lepage equivalent of a non standard problem. We will apply these 
considerations to our problem. The important thing to note is that, if a vari- 
ational problem has a covariant and contravariant Lepagean equivalent problem, 
then the latter can be considered as a kind of Hamilton- Pontryaguin's principle for 
the given variational problem; in fact, it is shown below that the canonical Lepagean 
equivalent problem associated to the variational problem underlying the Hamilton's 
principle gives rise to the classical Hamilton-Pontryaguin's principle (see example 
7). This will be our starting point for assigning a multisymplectic space to the 
variational problem we are dealing with. The bivariance ensures us that every ex- 
tremal has been taken into account in the new setting. Otherwise, namely, for non 
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contravariant canonical Lepage equivalent problems, some extremals for the original 
variational problem could be lost in the process. 

So let us suppose that we have the non standard problem defined by the following 
data 

'A ->• A x ->• M, 

x c Q* (A) , 
S[<r] :=J M (prcr)*(A). 

If I in this non standard problem has the required regularity (i.e. the "constant 
rank" hypothesis), then the variational problem (A — > M,T,X) will have a canoni- 
cal Lepage equivalent problem (\V X -> M, {0} , e) , and we can apply the scheme 

described above. In order to carry out this task locally, let us suppose as above that 
the fibers of the bundle / on an open set U C A can be written as 









7 >" 


• : <*ki ly' a l l 7 >" 


• » °4 1 



•••XI 



. • • • > < 



7 £ 17 



so that the prolongation structure I is (differentially) generated by 
X = ia\ ' : 1 < i < p, 1 < j < h 



on [/, where ( ol\, 



diff 



ICj) = : Xf]Q. j (U) and G Z x (A) for all i, j. Then 



J / diff 

if fTo := 7T (L7) c M, we can define for each 1 < I < p the numbers mi := dim (M) — l 
and the n-form on A v := U x Uo 0f =1 [A™ 1 (T*J7 )] efci will reads 



i=i \j=i 

where (/3, 1 „ 1 , • • • , ^ , • • • , ^ , • • • , /?m p p ) denotes sections of f\ mi (T*U ); the sub- 
script in these sections thus indicates their degree (in the exterior algebra sense). 
It will be shown in the examples below that, in many important cases, the Euler- 
Lagrange equations associated to the non standard problem defined by the data 

Au ^> Au -> AT, 
OcSl* (Au) , 

:= A/ (P rcr )* A ' 

has a family of solutions which is isomorphic to the family of solutions of the previous 
system; this means that in these cases the canonical Lepage equivalent problem is 
also contravariant. By proposition 1, the Euler-Lagrange eqs for an stationary 
section a G L ( Au ) are 



a* ( VjdX ) = 0, WeL 



(2) 



because there are no conditions for admisibility of variations; these sections are then 
integral sections for the EDS 



V"jdA : V G T VA 



L(7 



diff 



(3) 



We call this EDS the (local version of) Hamilton- Cartan EDS. 
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2.2.4. Some examples. We will illustrate with some examples how the standard 
approach fits into the scheme developed above, and moreover we will use it in 
variational problems with non standard description. 

Example 7. Let us start with the standard variational problem of mechanics. 
This scheme can be found in [18, 19]. Our Lepage equivalent problem will have as 
underlying bundle 

f\ (T* (M x TQ)) = T* (Rx TQ) = Ixlx T*TQ. 

By fixing U' C Q a coordinate neighborhood, we have that on f := ffi x TU' the 
EDS I is generated by the collection of forms 9 l := dq 1 — q l dt\ therefore the bundle 
JcT* (Ex TQ) will have the fibers 

Thus I 1 = K 1 (the elements in I are of first degree in the generators); moreover 

T* (K x TQ) c Z 1 {Rx TQ) 
implying K 1 = K 1 n Z x (R x TQ) = I 1 . In terms of our local description we obtain 



that W } 



U = U x R n , where 



(t,q,q) 



Ldt + a J 0j : Oi 6 R, * = 1, • • • , n \ V (t, q\ «*) G 17. 



The o's in the previous formula have an important meaning concerning the manifold 
A:=Ix (TQ © T*Q) with the 1-form A e il 1 fA) defined through A := 

V / (t,v.a) 

— [L (a) — a (v)] dt + <df (where <df is the canonical 1-form on T*Q), as the 

a 

following lemma shows. 

Lemma 2.9. The map Pl : W\ — > A given locally as 

(n \ / n r, n \ 

is a diffeomorphism, and P£A = f2. 

We can work then on the simpler space A = R x (TQ T*Q), where locally 



X = J2pi (dq l - q l dt) - Ldt. 



i=i 



The canonical Lepage equivalent problem consist in finding sections a of R x 
(TQ © T*Q) which extremize 

S [a] = j a* (^JPi (dq 1 - q l dt) - Ldt^j ; 
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the equations for an stationary section a : t i-> (q (t) , q (t) ,p (t)) are then 

'dq* - tfdt = 
dpi - L qi dt = 
( Pi - L#) dt = 0. 

We identify here the Hamilton-Pontryaguin variational principle. 

Example 8. The standard first order field theories deals with sections a of a bundle 
ir : F — > M with typical fiber V; a lagrangian density L G C°° (J 1 (7r)J allow us to 
define the action 



S [a] := I (Lo prcr) dx° A ■ ■ ■ A d; 

J M 

= I (L o jV) da; A • • • A dx 



,71-1 



ra-1 



according to considerations previously made. Let U C F be an adapted coordinate 
chart with coordinates (x fc , u^) ; then X on U := J 1 (n) | W has the generators 

{6> Q := du a - u%dx k : a = 1, • • ■ , dim^} C ft 1 (u 



where the Einstein's summation convention was asummed. So in this case we will 
have 



W 



U 



{a + fc- 1 a e a ■. 6 z'o - 1 (u) } 



by using the F\ map it can be written 



u~u< 



3dimV 



(where Uq := tt (U)) and this is the local description that we will adopt. Therefore, 
we can consider the Cartan form defined on A := U ® Z with local coordinates 



©dimV 



denotes the set of diml/ forms of 



(x k ;u a ,u%,m a ); here Z := f\ n 1 (T*U ) 
degree n — 1 on M. So we can write 

A := m a A (du a - u%dx k ) - Ldx° A ■ ■ ■ A dx 12 ' 1 

for the Cartan form. If the variations are (0; Su a , Su^, Sm a ) then the Hamilton- 
Cartan equations reads 



(0;5u a ,5u%,5m a )jd\ = 



dL 



6u a 



dL 



5u? I dx° A ■ ■ ■ A dx 



n-l 



A (du a - u%dx k ) 



Su a dr 



m a 5u%dx 



Now we can define the collection of (n — l)-forms (r a ) on M via 



3T 

— dx° A ■ ■ ■ A dx 71 ' 1 = T a A da 



Va, k, 



from which we obtain that the solutions 

x k i->- (x k ; u a (ar) , u% (x) , m a (x)) 
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must satisfies 

r a — m a = 
du a - ul&x k = 

dm a - ( J^) dx° A • • • A dx"- 1 = 0. 
This system projects onto the local version of Euler-Lagrange equations for first 
order field theory. The space of (multi)momenta is locally U © Z. 

Example 9 (Electromagnetism - Cont.). The first true hamiltonian non standard 
theory which we want to deal with is electromagnetism; in this case we have that 
A := A 2 (T*M)©A and A = p\ (*9 2 ) A (dp^Gi -p%& 2 ) + \, beings the projection 
in the i-th summand of A and A the pullback of the corresponding form on A. 
As before, let us suppose that we indicate by (P, F, A) 6 A an element of this 
(multi)phase space; then a section a : x M- (P (x) ,F(x),A (x)) will be a solution 
for the eqs 2 iff 

' F = dA 
P = *F 
dP = 0. 

Thus the space of momenta coincides with the space of velocities. 

Example 10 (Poisson sigma models). The Poisson sigma model [G, 3, 4, 26] can 
be analysed from this point of view. This is an interesting system, because its usual 
description (cf. the references above) is non standard, in the sense that it deals 
with variations without a prolongation structure. Concretely, for a surface E and a 
Poisson manifold (M, it) we have the fibration 

2 

A := f\ (T* (M x £)) -» S x M -> S. 

Then the open set S := Z\ (E x M ) — Z% (E x M) is a subbundle; let us call 
v : S — > E the restriction of the above projection. The non standard variational 



problem equivalent to the Poisson sigma model is the triple \ S — > E, 0, A J , where 

in the local coordinates (£ Q , x^, rjp u ) associated to a coordinate system on M x E 
via 

a G 5| (a . )£) i-> -ripvdZ A dx v , 

the 2-form A reads 

1 

(Cx, v ) ' 2' 
The Poisson structure on M defines a bundle map II : T*M — !> TM and there- 
fore induces another bundle map (which we denote with the same symbol) II : 
T* (M x E) — > TM x T*E; then the form A can be written in global terms as 



A 



- mv dx v A d^ + -Tr^V^dr A d ^ 



A 



where {•, •) is the pairing of an element of TM x T*E with an element of T* (M x E) ~ 
T*M x T*E and <d 2 is the canonical 2-form (restricted to S). The local generators 
of the Hamilton-Cartan EDS are in this case 

fd^-Tr^^de*, 

\d Vaa a de* - i (d^) Va ^d^ a d^. 

Thus these systems are naturally formulated in a multisymplectic space. 
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3. Slicings and presymplectic structures. Given a Lepagean variational prob- 
lem (that is, a non standard problem with trivial prolongation structure) and a 
decomposition of the underlying bundle in a family of hypersurfaces (i.e., the con- 
stant time slices), the set of sections of a fixed slice can be dressed up with a 
presymplectic structure and a hamiltonian, under convenient regularity hypothe- 
sis; these data yields through Gotay-Nester algorithm to a set of constraints on 
the sections under consideration. The procedure is explained with some extent in 
the following paragraphs; additionally, it is settled the relationship between the ex- 
tremals of the variational problem and the solutions of the Hamilton eqs associated 
to the hamiltonian and the presymplectic structure on the space of sections of the 
fixed slice. Finally, it is shown how the usual Dirac theory of constraints for a field 
theory can be derived in this setting. 

3.1. Preliminaries. The purpose now is to set the framework allowing us to relate 
the solutions of the Hamilton eqs on n — 1-sections with regular and ordinary n — 1- 
integral elements of some associated EDS. The initial data is a bundle F —> B and 
a n-form A (n — dim_B) on F. Then the Euler-Lagrange eqs for the action 



J B 

characterizes the sections a : B — >• F which extremizes it. These sections are integral 
for the EDS (we assume that dB = 0) 



We would like to introduce a presymplectic manifold with a hamiltonian, whose 
solution curves are in one-to one correspondence with the integral sections of I. In 
this vein it is necessary to introduce the concept of slicings. 

3.2. Slicings. The first thing that we need [16] is a compatible slicing (sf,sb) of 
the bundle of fields: It consist of a pair of diffeomorphisms 

s F : ffi x K -> F 
s B : K x S -)• B 
making commutative the diagram 

E x K — F —- F 



Note 1. Although not strictly necessary, we will consider that the factor £ in the 
decomposition of the base manifold is compact. This makes sense to definitions like 
10 and 11 below, and turns the 2- form automatically closed. 

Let us define F T := sp ({t} x K) and S r := sb {{t} x E) for each r £ R. A 
useful consequence of this property is that the space of sections V (F) admits a 
trivialization; in fact, defining for each rfl the embedding 




I:= (ZjdA : Z e T (VF)) 



diff • 



(4) 



IxE 



B 



i T : S — > B : x n> sb (t, x) 



we have that the map 



s T{F) : K x r (A") -> T (F) : (r, a K ) i— > s F o a K o (i T ) 



-l 
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makes it works (it is well-defined because of the compatibility of the slicing, meaning 
in particular that Im (sf \ S t ) C A r .) 

3.3. (Pre)symplectic structure on the space T(K). For each r G R we have 
the diffeomorphism i T : A — ¥ F T : x t— > sp (t, x) so we can define [10, 8] the following 
r-family of forms on A 

n T -. = i* T (d\\ f t ) e n n+1 (K) 

H T : = i*(d jX\F T ) G ft"" 1 (A) 

where the symbol do is used in order to denote the canonical vector field associated 
to the r-direction in both spaces. These forms can be used to define the following 
data on T (K) for each r G R: The 2-form uj t G fl 2 (T (A)) and the function 
H T G C*°° (r (A)), defined through 

u T \ 7 (Vi,V 2 ) : = ^7*(^r (^1,^2)) 

Hr{l). = fl*{Ur). 

Here 7 G T (A) and Vi,V 2 G Y{^*(VF)) = T 7 r (A); as before, t> G r(M) 
denotes an extension to a vector field on A" for the element V G T (7* (VA")) . With 
this data we can define a Hamiltonian system on T (A). 

3.4. Forms on T (A) and forms on K. We want to establish a useful result 
concerning forms on a space of sections; that is, by taking on a bundle F B an 
n + /c-form <f> G Q, n+k (A) , fc G Z, n = dimS, we can define the following fc-form i> 
on T (A) via 

$| CT (Fi, ■ • ■ , 14) := j o* (Vu ■ ■ ■ V k j4>) ■ (6) 

-n + k 

Definition 3.1. For each k G Z define the bundle Z£ +fc (A) ^ A such that the 
fiber 011 e £ f is given by 

{n+k ~| 
a G /\ (A* A) : l/i j ■ • • Ffc+ua = Vft, • • • , Vfc+i G T4A > . 

This set can be called the space of k-semibasic n + k-forms. We will say that a 
nonzero form a belongs to Z* properly iff a G Z* but a ^ (A). 

Then we can prove the next result. 

Proposition 3. Let us suppose that (j) G T (Z£ +l (A)) properly (I < k), and let $ 
be as in Eq. 6. Then $ = iff <f> = 0. 

Note 2. The previous proposition can be rephrased as follows: The map <fi 1— > $ is 
injective when restricted to T (A)). 

Proof. It is clear that </> = implies <f> = 0. For the converse implication, $ = 
implies that j B a* (Vij ■ ■ -Vijifij = for all a G T (A) and Vi, • • ■ , V/ vertical 

vectors defined on Imcr; therefore V\ j • • • Vi -i(f) = for all collection of I vertical 
vectors, and it means that cf> G (A). Therefore <f> = 0. □ 
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For each relwe can define the following EDS on K 

1 T : = i* T (l\F r ), (7) 

a fundamental object in the subsequent sections, because it is closely related to the 
Dirac constraints of the theory. 

Finally we state a very important result, because it relates the solutions (if any) of 
the presymplectic dynamical system with the extremals of the underlying variational 
problem. 

Proposition 4. Let X be a n-form on F such that X € Z™ (F) properly. Then for 
each Tfl, a pair (cr T , X T ) G T aT T (K) satisfies the Hamilton equations for sections 

(XrJWr)U = dffrU (8) 

iff it satifies the following equation 

o-; (v_i (x t m t - <m T )) = 0, Wer « (vk)) . (9) 

Proof. The key is to realize that the n-form X T jf2 r — d"H T belongs to Z^ (F) . On 
the other side it was shown in proposition 3 that there are no non zero n-forms fj 
in Z% (F) such that the map 7 4 J s a* (V j/3) is zero. □ 

The meaning of this proposition must be clarified: It says that the pair (<r T , X T ) 
is a solution of the Hamilton equations on T (K) defined by u T and H T if and only if 
V T := da + X T belongs to the polar space of the n — 1-integral element ov, (TS) for 
the Hamilton-Cartan EDS in this context, Eq. 4. Then if there exists a regular flag 
for this EDS with last terms Oc ■•■ C d T , (TS) C (<t t * (TE) , V T ), the n-integral 
manifold passing through it will generate solutions for the Hamilton equations, and 
conversely. The following sections will use this remarkable relation. 



3.5. Presymplectic structures on a multimomentum space. Let us now ap- 
ply the setting of the previous section to our multisymplectic space. In particular, 
let us suppose that there exists a compatible slicing of the bundle A — > M: So be- 
sides of the spacetime decomposition M = Ex E (as above, it will be assumed that E 
is a compact manifold), we have a decomposition A = RxI such that A T = {r} x L 
for every t£1; thus p : L — > E is a fibration that makes commutative the following 
diagram 



x L ~ A 



idxp 



x E ~ M 



P 2 



L 



In such a case (cf. the previous section) we can define a presymplectic structure on 
the space of sections T (iA together with a function on it such that the solutions of 
the classical mechanics system defined by these data induce solutions for the EDS 
generated by 3. That is, if we define 



L -> A r : / h-> (t, I) 
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c^A |A T 



we can use these structures to define the forms 

t(l 



dX 



so, we have H T <S 
dependence) defined by 



and 



e n a r £ 



H T (a) 



<?* {Hr) , 



(10) 

(11) 

(which preserves a r- 

(12) 
(13) 



for every pair X T , Y T of vertical vector fields that cover a, i.e., elements of T a T 

In this setting the proposition 4 is true, and we have to our disposal (under adequate 
regularity conditions, i.e., if the foliation induced by the slicing is regular in the sense 
of definition 4.3 below) the identification between the dynamical equations on the 
space of sections and the Hamilton-Cartan EDS. 

Now it will be necessary to introduce further simplifications; in particular, we take 
Cg X = 0. Using the slicing introduced on A and the Hamilton-Cartan EDS, given 
by 3, we can define a parameterized EDS I T according to the Eq. 7. Under this 
simplification it is true the following result. 



Lemma 3.2. If Cq X = 0, then H T = H T * =: H andco T 
Moreover, the EDS X T does not depends on t. 



uj for all t, t' € 



In particular, this lemma allow us to choose a slice to work with. 
Finally we want to show how the usual theory of constraints in a field theory can be 
obtained from this scheme. In the following example we set up the relation between 
the presymplectic structures defined above and the symplectic structures defined in 
the usual approach to (first order) field theories. 

Example 11 (First order field theories - Cont). Our aim is to calculate the dy- 
namical system associated to the multisymplectic structure for field theory found 
in example 8; locally we have that A = m a A (du Q — u^dx k ) — Ldx° A • • ■ A da;™ -1 , 



dA = dm a A (du a - u%dx k ) — dL A dx° A • • • A dx" -1 - (-if ^ m a A du% A dx k , 
and taking into account that m a G il™ -1 (M), we will have that 



«dm a A dx k )\ No = (m a A du% A dx k ) \ Nq = 

where N := (p e U ® Z : x° (P) = oj. Then if N := {m E M : x° (to) = 0}, by 
using Eq. 13, the presymplectic form on T (No) will be 



uj x \ a (X, Y) = (Sm°Jv a - Sn°Ju a ) dx 1 A • • • A da- 



Tl-l 



N 



and X := (0; Su a , 611%, 5m a ) , Y := (0; 5v a , Sv%, <5n a ) indicates a pair of tangent 
vectors at a S r(iVo)- By introducing the coordinates 

m a : = m k a d k x 7 
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with the symbol d^x defined through 

d k x := d x k j (dx° A • • • A dx™" 1 ) = (-if dx° A • • • A dx^ A • • • A da 
then we will see that for the restriction of the n — 1-form 

d x oA = (d x ojm a ) A (du a - <dx fe ) - Ldx 1 A • • • A da;"" 1 - 
to No, it can be proved that 



d T o j A 



= (d x ojm a )\ No A (du a - u a A dx A ) — [L - (-1)™ ti a m°J dx 1 A • • • A dx™" 1 , 
where the capital letter indices takes values in the range 1, • • ■ , n — 1. Finally 
(d x o jra„) A dx B = (c^ojd^x) A dx B 

= (-lf^mfdx 1 A-.-Adx 11 " 1 
and from Eq. 12, the hamiltonian function H x £ C°° (T (iVo)) will becomes 



H* (a) 



a* cLojA 



iVn 



{(-l)"" 1 m B [(d B u a ) -u a B ]-L- (-1)"" 1 } dx 1 A • • • A dx"" 1 . 



We are now ready to make contact with the usual theory of constraints. In order 
to do that, it is necessary to calculate the primary constraints associated to the 
dynamical data found above. This can be achieved by realizing that keru;^ = 
{ (0; 0, Su^, Sm A ) }, so the primary constraints for the dynamical system will be 

= diP {5ut) 

dL 



N 



(-!)"< 



Su^dx 1 A ■•• Ada;™" 1 , 



= dH x (Sm£) 
= [ (-1)" -1 [(d B u a ) - u B ]8m B dx 1 A- • ■ A dx 

JNn 



n-1 



'N 

Let us now define the set of primary constraints 



Ci 



-1 m 



a = 1, • • • dim]/, k = 0, • • • , n — 1, B = 1, • • • , ?i — 1 



then Mi C T (No), the zero locus of C\, is a presymplectic manifold with hamilton- 
ian, both structures being defined through restriction. Let No be the trivial bundle 
iV x M m x M" 1 '"- 1 ' x W n x R m on N with coordinates (x; Q , tt q ); then 

the map II : iVo — > Afo defined through 



n : < 



' (j) a 


= U a , 


<t>% 










= (-1 
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has the following properties. 

Lemma 3.3. Let L\ € C°° (A/"o) be the unique map such that L± o II = L. The set 
Mi := II o (Mi) is composed by the sections of ATq given by 

a : x (x; <j> a (x) , (j>% (x) , <j% (x) , ir a (x)) 

such that 

din. , a ^ a 



Ox^ 

If Mi is a manifold, the the presymplectic form defined through 



/ dn a A dcj) a A dx 1 A • ■ • A dx n 

JNn 



I N a 

will verify that (IIo) uii = uj x \ Mi ; also the hamiltonian function 

Hi (a) = [ a* Oo7T Q - Li) dx 1 A • • • A dx"" 1 

will satisfy (IIo) Hi = H x \ Mi . Due to the injectivity o/(IIo)* ; the structures are 
uniquely determined by these requeriments. 

Then the dynamical problem reduces to solve 

X Hi auji = AH i 

on each o \ x i — y (xj (j) a (x) , <\)°^ (x) , 0g (x) , 7r Q (x)) such that 

_ dLi _ dr 

and this is the basic scheme in the usual theory of constraints for a (first order) 
field theory; therefore, the primary constraints in the standard approach are related 
through II with the primary constraints found in the non standard viewpoint. 

4. Constraints and EDS. The aim in this section is to prove the main result of 
the article: The fact that, under certain regularity hypothesis (cf. definition 4.3 
below), it is possible to describe the set of constraints arising from the application 
of the Gotay-Nester algorithm on the dynamical data 12, 13 as a set of generators of 
the EDS induced on the chosen slice by the (prolongation of the) Hamilton-Cartan 
EDS. 

4.1. Useful characterization for constraints. As we said above, we want to 
relate the Dirac constraints arising from a field theory (in the non standard setting) 
with a EDS associated to its eqs of motion. In order to achieve this, it is necessary to 
adopt a meaningful picture for these constraints. So we use the following description 
for the constraints obtained through the Gotay-Nester algorithm [16]. 

Definition 4.1. Given a presymplectic manifold (N,u) with a hamiltonian H G 
C°° (N), the constraint submanifold C C N is the maximal submanifold of N 
according to property 

dH (TC^) = (14) 

where TC 1 ^ denotes the (perhaps singular) subbundle of TN composed by the 
symplectic complements of the fibres of the subbundle TC . A submanifold in M 
satisfying property 14 is said to be invariant respect to the dynamics defined by 
(lj,H). 
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Because we need to know the tangent spaces for submanifolds T (/C) C T ^Z/j 
composed of integral sections for some EDS /C, we prove the next lemma. 

Lemma 4.2. If the set T (JC) , with K, C fi* (^Lj an EDS, is a true submanifold in 

r (^Lj , then its tangent space in a £ T (JC) is the vector space 

T a T ()C) = {ler (a* (VL\ ) : a* (C t K) = o} , 

where the hat denotes some extension to X[L) of an element of T (a* I VL 



Proof. Under our hypothesis, we just need to derive along a curve to obtain the 
lemma. □ 

4.2. Admissible sections and involution. We start with a (non standard) vari- 
ational problem on a bundle A — > M, where dimM = n; by assuming that there 
exists a (bivariant) Lepage-equivalent problem, we pass to a non restricted varia- 
tional problem on a new bundle A — » M. The variational equations defines here an 
EDS Icfi' (j^J ; the integral sections for this EDS are the classical solutions to 
our field theory. So we have the EDS 

(X, n := dx 1 A • • • A dx n ) 
generated by the Hamilton-Cartan equations, i.e. 

X := (VjdX :VgT (vAj ^ , (15) 

and by using the Cartan-Kuranishi theorem we found a bundle A' — > M, a submer- 
sion II : A' — > A and an EDS (Z',S1') C O* ^A'^ which is n-involutive and verifies 
that 

il (v n (i',n')) = v n (i,n). 

Remark 1. It is important to pointing out the following facts. 

• The original EDS does not work if there exists n — 1-integral manifolds on 
each r-slice whose possible thickening directions rests in the same slice. The 
way to avoid this annoying fact is to select in V n (I) only those integral el- 
ements satisfying the independence condition; this subset of sections can be 
described as n- integral submanifolds for an involutive EDS X 1 , obtained from 
X by performing enough prolongations. 

• It is assumed that there are no 0- forms to take care of in X and X'\ if it occurs, 
one must to include them into the bundles A and A', that is, redefine these 
sets by taking into account that the new 0-forms annihilates on them. The 
corresponding EDSs are the pullback of the original ones. 

• We know that A — > M admits a compatible slicing; we will suppose further 
that the same is true for A' — > M with A' ~ R x I/. This diffeomorphism will 
be denoted as s^, . Moreover, it will be assumed the existence of a submersion 
IIo : L' — > L such that the following diagram is commutative 

n 

A' A 



V L 
iio 
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• For each r € M it is tempting to define the EDS 
X' T :=i*(x>\A> T ) CO* (£') 

where A' T := i T ^{r} x L'J and 

i T : L' <-+ A' : Z h> s^, (t, /) , 

claiming that 

r n _ 1 (i' T ,uj') = i*(r n (i',n')) (ie) 

where u/ := <9 jO'. It is obvious that r„_i {T T ,u') 3 i* (r„ {X',0,')). The 
n-involutivity of X' ensures that every n-integral element is tangent to some 
solution; however, the opposite inclusion r„_i(Z^,u/) C i* (T n (X', f2')) is 
achieved if every n— 1-integral element of X' T can be thickened out (in a regular 
fashion) to a n-integral element for I'. The requeriment of involutivity is not 
enough; instead, we need the regularity of every n — 1-integral element of X' T . 

Note 3 (System where 16 is not verified). There exists an almost canonical coun- 
terexample to this claim: The EDS describing orthogonal coordinates in three di- 
mensions. See example 3.2 and theorem 3.3 in [">]. 

Due to the previous discussion, we need to introduce a new definition. 

Definition 4.3. Let X be an involutive EDS on M and let us suppose that M ~ 
R x N for some manifold N; as above, for each r 6 1 let us define the following 
EDS on N 

T T := i* (1\ {t} x N) . 

We say that the induced foliation is regular with respect to the EDS I if every n — 1- 
integral element of Xr, r G K, is regular. A slicing inducing a regular folitation will 
be called regular. 

The regularity notion for a foliation ensures that one can extend any n — 1-integral 
manifold for the EDS contained in a leaf to a n-integral manifold; for suitable 
independence conditions the required involutivity ensures that the extension is not 
included in that leaf. So we can state the following result. 

Proposition 5. If I is an involutive EDS on P, and sp : P — > Q xl is a regular 
slicing for I, then 

r„_ 1 (z r> w) = i; (r„(z,n)), 

where lo := dojfl and 8q is the vector field pointing in the ^-direction of the slicing. 
Proof. The hard inclusion is 

r„_i(i;,w) c»; (r n (i,n)). 

Because of the regularity, it follows from the Cartan-Kahler theorem that any n — 1- 
integral section for I T can be extended to a n-integral section for X. □ 

It remains to show some revelant properties of the solutions of an EDS. In order 
to do that, it is important to establish the following fact. 

Lemma 4.4. Let X be an EDS on the manifold M , N C M an integral submanifold 
and leT (T N M) such that 

(Xja)\N = 
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for all a el. If X e X (M) is any extension for X , then 

Proof. If i : N '—t M is the canonical injection, we have that 

i* (Xjcx) = i* (Xja) , 



so 

(C^a)\N = i* (x_ida + d(Xjaj^ 

= i* (Ajda) + di* (Xja) 
= 

for all a el. □ 

By using the definition 7 we obtain an EDS (Z^.,u/ := do-ift') on L'. As we will 
see in the next sections, this exterior system gives another description for the Dirac 
constraints, namely, as (algebraic) generators of an EDS. 

4.3. The constraint submanifold as an space of sections. 

4.3.1. An outline. Here we describe the essential elements involved in proving that 
the constraint submanifold coincides with the set of maps iloor„_i (X' r , uj'). Because 
this proof uses the characterization of C as the maximal invariant submanifold of 
the set of section of the bundle L, it is necessary to prove that 

1. Bio o r„_i (1' t ,uj') is invariant, and 

2. <7cn oIVi(X;,a/). 

It seems that the first condition is consequence of the fact that the original PDEs 
9 are satisfied by the elements of this set. The second condition can be verified by 
building a n-section from an— 1-section in C via the formula 

a : (t,x) (t,oy (a;)) 
realizing that it is an integral section for (Z',f2'), and then restricting to a r-slice. 
It was previously stated that a submanifold Q C P in our phase space T (j^j is 
invariant iff 

dH {TQ^) = 0. 

According to a result of Marsden et al [16], the invariance property is equivalent to 
the existence of tangent solutions to the Hamilton eqs restricted to Q; that is, for 
every q e Q the system 

{X H ^)\ q = dH\ q 

has a solution Xjj\ q belonging to T q Q. In the first part of the proof the construction 
is performed in order to ensure that the submanifold defined by the EDS admits 
tangent solutions to the Hamilton eqs, implying the invariance. In the second part 
the strategy is to prove that our submanifold is maximal with respect to the property 
of invariance; now the characterization of invariance is used in the opposite direction, 
that is, it is initially supposed that certain submanifold is invariant, and by using the 
characterization, arriving to the conclusion that the Hamilton eqs. admits tangent 
solutions to this submanifold. Finally it can be related to the manifold defined by 
the EDS because of the interpretation of the Hamilton eqs as generators of the EDS. 
After this warm-up, we are then ready to formulate the theorem with its proof. 
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4.3.2. The main theorem and its proof. There are several elements to be included 
in the final proof of the main theorem of the paper. In the previous paragraphs 
it was argued the necessity of having a submersion IT : A' —> A between bundles 
fitting in the diagram 

n 

A' - A 



M 

Next it is assumed that there exists a compatible (with respect to the both bundle 
structures) slicing of M which is well-behaved with the maps in this diagram. That 
is, there exists a map EL} : L' — Y L making commutative the following 



A' 



V 



n xid 



A 



L x 



The slicing on A' allows us to define the EDS 

K ■= % (l\ L' x {r}) , 

which is nothing but the original involutive EDS X' restricted to the leaf V x {r}. 
In order to ensures that V n -\ (X' T ) C V^_i (X'), it is required furthermore that the 
slicing be regular with respect to the EDS I'. Now we define 

v{x' T ) ^riooiVi^V), 

and the main theorem can be stated as follows. 

Theorem 4.5. Let us suppose that the EDS X defined through 15 is r-invariant, 
namely 

C do X c X. 

Then the constraint submanifold C in the presymplectic space of sections is equal to 
T(X>). 

Proof. The constraint submanifold C can be characterized as the maximal invariant 
submanifold of the set of sections of L [16], where invariance for J C T (j^j means 
that 

AH (TJ X ) = 0, 

being H the Hamiltonian function for our system and the orthogonal complement 
is taken with respect to the presymplectic structure in the space of sections. It 
can be proved (prop. 6.9 in [16]) that the invariance of a submanifold J in a 
presymplectic manifold is equivalent to the fact that through any point of J there 
exists solutions to the Hamilton equations tangent to this submanifold. So let us 
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take cr T G r (X' T ); according to the very definition of this space, there exists a unique 
< G r„_i (I' T ,u') such that 

cr r = n o a' T . 

Because of the regularity of the submanifold Im {a' T ) C A', by Cartan-Kahler 
theorem there exists a (non necessarily unique!) integral section a' G T n (I', f2') 

such that cr'l ^1/ x {r}^ = cr^.. If 9o denotes alternatively the vector fields defined 

through 



0, 



0|(A,r) 



(m,r) 



= (*a<); 

= ( S A/), 



(A,r) 



(m,r) V dr 



then we can define the section X' T G T ( (cr' T )* \ VL' 

(a% (d Q ) = 8 Q + i„ (X' T ) 



on A' 
on M 

t(l' 



(17) 



Thus the expression X T := IIo* °X' T is a well defined section of the pullback bundle 

(ov)* (vlJ, because the prolongation procedure and the proposition 5 ensures us 

that the map IIo ° (•) between n — 1-sections of the relevant bundles is one-to-one. 
Furthermore, by projecting both sides of 17 along II* o (•) we obtain that 



er* (do) = d + i T * (X T ) 



(18) 



where a := II o a' is a n-integral section for X on A, due to prolongation property. 
In particular, ct* |/ t s (T m M do) C TA is a n-integral element for I contain- 
ing the subspace cr T *| m (T m M), which is n — 1-integral for X T . Then ct* (do) G 
H (cr T *\m (T m M)), so we must have that 



(d )ja)\(a T *\ m (T m M)) 
By using eq. 15, we obtain that X T satisfies 

a* (Vj (X t M - cm)) = 0, 



0, 



Va G 1. 



yv g r ( cr* Ivl 



meaning through proposition 4 that X T is solution of the Hamilton eqs. associated 
to H. It remains to show that this solution is tangent to T (T' T ); according to lemma 
4.2, it is equivalent to show that 







for some extension X T G X (Lj of X T . Because of lemma 4.4 and taking into 
account that o> is a n — 1-dimensional integral submanifold of I, we conclude that 



K (£ x i) = o, 



(19) 



where X := do + i T yX T J, and X T G X (Lj is some extension to X T . The r- 
invariance 

C dn l G 1 
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implies that 2 

<(A„(a t ) I )= ' ( 2 °) 

and furthermore, because X T (to) £ VL for all m £ M, for each a £ I we will have 
that 



a* T {i„(X T )^a) = a* [i T * [X T )j [a\ L X {r}j 
= cr* (x t j (i* a| L x {r }H . 
Then the Cartan's Magic Formula allows us to rewrite the expression 20 as 

< UxZt) = 0, 



which is the desired tangency condition. Therefore we have shown that T (I' T ) is an 
invariant submanifold of T (i/j . 

Our next task is to show that T (I' T ) is maximal among all the invariant submanifolds 
of the space of sections T ^Lj . So let us take Q C T (^Lj an invariant submanifold. 

Then for every <jq £ Q there exists Xq £ T ^CTq (vLj J such that 

• X £ T ao Q, and 

• a* (Vj (X M - dH)) = for all V £ T (a* (vLj ) . 

Assuming that any vector field Z £ X (Q) assigning to every section ao a vector X 
with these characteristics has integral curves 3 r > a T passing through <tq £ Q when 
r = 0, we are able to build the n-section a : (r, x) i— > (r, a T (x)) of A. This section is 
n-integral for I because of the Hamilton equations. In fact, both conditions above 
implies that er r * (T X T,)® (<9o| T ) (r belonging to some neighborhood of 0) annihilates 
the set of forms 

S := {VjdA : V £ T (vAH ; 

because S C 1, we have that (<S) diff = I (by the definition of I, see Hamilton-Cartan 
EDS, equation 3) and thus V n ((<S) difi ) = V n (I). Then by using the prolongation 
properties we can find an integral section a' which is n-integral for I' and Ilocr' = a\ 
moreover, 

i o(T = nooTa-'IAoJ 

and a\ A' Q £ r n _i (Xq), and then a £ n o r n _i (2J) = V (I' ). □ 

In the following section we will use the relation settled by this theorem between 
Dirac constraints and the EDS X T in some interesting examples. 

5. Examples. This section contains examples where the techniques developed in 
the work will be applied. The first two of them deals with variational problems re- 
lated to field theories; the last two shows some applications with more mathematical 
taste. 



2 Abuse of language: The symbol i r * (^tJ in the following equation represents the vector field 
(/, s) i— > i s * (^X s (m)^ on A ~ L X R, where the number r refers to the second factor in this 
decomposition. 

3 We hide behind this assumption some issues concerning hard analysis. 
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5.1. Techniques in the EDS treatment. Before to start with, it is necessary to 
mention some specific issues related to the verification of involutivity for a system 
with independence condition. Namely, it consists of two parts: 

1. The verification that the reduced characters are equal to the original charac- 
ters. 

2. The verification of the equality hypothesis of the second part of the Cartan's 
test. 

If any of these fails, we use some of the tools discussed in the appendix B in order to 
find an EDS with the same integral manifolds (in some dimension) which make them 
hold. In order to work with these conditions, let us suppose that E = (vi , • • • , v„) C 
T X X is an n-integral element of an EDS (I, U3\ A • • • A u> n ) on a manifold X and we 
take the flag induced by the subspaces 

E : = 0, 

E k ■ = (vi,---,v fe ), k = 1, • • • ,n - 1; 

for each k, the linear system 23 gives us the space H (Ek) as the kernel of a linear 
operator M (vi, • • • , V&) : T X X — > T*X. Then we have that the Cartan characters 4 
are 

c fc = rank M (vi, • • ■ , v fe ) . 

For the calculation of the reduced Cartan characters, we must write the operators 
M (vi, • • • , Vfe) in terms of a basis {lo\, ■ ■ - ,u n , • • • } C T*X and its dual: The 
reduced Cartan characters Ck are the ranks of the matrices obtained through this 
procedure, after deleting the columns corresponding to the elements u>\, • • • ,u„ of 
the chosen basis. 

5.2. Electromagnetism. We apply our results in the case of the classical theory 
of electromagnetic field; as we know, the EDS associated to their Hamilton-Cartan 
equations is defined on A = A 2 (T*M)®/\ 2 (T*M)®T*M with coordinates (P, F, A) 

by 

I := (F — dA, P - *F, dP, dF, dP - d (*F)) alg 

with the independence condition Vl := (jp% o tm)* (wm) 7^ 0, where ojm £ Q A (M) 
is some volume form defined on the space-time (by supposing M orientable). For 
this system the reduced characters £2 and C3 are less than the corresponding Cartan 
characters; this is because of the appeareance of the form P — *F in the system, 
which reduces to a collection of functions on A when it is evaluated on integral 
elements satisfying the independence condition. In order to circumvect this, we will 
deal with the associated EDS on A := A 2 (T*M) © T*M defined through 

J:= <F-dA,d(*F),dF) alg ; 

it is the EDS induced by pullback of our original one to the submanifold described by 
the equation P = *F. Anyway, by using the identification of the vertical directions 
with the fibres in our bundle, we use the following vectors 

V» := (d^-F^A^) 

as the basis inducing the flag 

E := {0} C Pi := (V 1 ) C E 2 := (V\ V 2 ) C P 3 := (V\V 2 , V 3 ) C E 



4 These are not the true Cartan characters as they are defined in the literature, although they 
are closely related. 
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in E := (V Q ,--- ,V 3 ) C T^^A whose Cartan characters are being calculated. 
Concretely 

• Because := J HQ 1 (A) = {0} we obtain that 

H (E ) = {« g T (F , A) A :vj(f> = for all e } 

= T(F,A)^ 

and so Co = codim H (Eq) = 0. 

• In this case 

H(Ei) = {ve T {FtA) A : vj (w-i<f>) = for all 4> E J {2) and w £ Ei\ , 

so u s H (Ei) if and only if 

vj (duF-A 1 -du&A) = 0; 

then a = rank A/ (V 1 ) =1. 

• Now 

H (E 2 ) = |v G T( F ,A)h- ■ UJ (wiJ(f>) = and uj (to2JW3J^') = 

for all (f> e J {2 \ip £ j' 3 ' and Wt, w 2 , w 3 € £2} 

and then v £ H (E2) iff it is solution for the system 

vj (<9uF - A 1 - dudA) = 

vj (d 2 jF ~A 2 - d 2 jdA) = 

vj(du(*F 2 ) -d 2 j(*F 1 ) +dud 2 jd(*F)) = 

vj (duF 2 - 8 2 jF 1 + dud 2 jdF) = 0; 

this means that c 2 = 4. 

• Accordingly, the equations determining H (E3 ) will be 



VJ 
Vj 







vj (duF -A 1 - dudA) 


= 






vj (8 2 jF -A 2 - d 2 jdA) 


= 






vj (d 3 jF -A 3 - d 3 jdA) 


= 


(du 


{*!*) 


1 -du^F 1 ) +dud 2 jd{*F)) 


= 


(du 




1 -dsj^F 1 ) +dud 3 jd(*F)) 


= 


(0 2 j 




1 -d 3 j(*F 2 ) +d 2 jd 3 jd(*F)) 


= 




VJ 


(duF 2 - d^F 1 + dud 2 jdF) 


= 




Vj 


(duF 3 - OsjF 1 + dud 3 jdF) 


= 




Vj 


(0 2 jF 3 - 8 3 jF 2 + d 2 jd 3 jdF) 


= 



the associated Cartan character being 

c 3 = 9 

because of the Hodge operator. 
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Let us now calculate the codimension of V4 (J) C G4 (TA); as before, it will be the 
rank of the linear operator behind the linear system 5 

V^jV^j (F - dA) = 

FWj(d(*F)) = 

V jV"jV a j (dF) = V/i, i/, a- = 0, • • ■ , 3. 

Using the symmetries of the contraction operator, we see that the a priori indepen- 
dent equations in this set are determined by the set of pairs of indices 

(0,1), (0,2), (0,3), (1,2), (1,3), (2, 3) 

and the set of triples 

(0,1, 2), (0,1, 3), (0,2, 3), (1,2, 3). 

The linear operator is constant, consequently it has maximal rank everywhere; so 
we have that codimsV^ (J) = # equations = 14. Thus the hypothesis of the 
Cartan's test are satisfied, and the EDS J will be involutive. By using the slicing 
:r° = constant we can identify M ~ K x L 

2 2 
K x o ~ (R x L) ® T*L ® /\ (T*L) © f\ (T*L) 

with coordinates (ao, a, e, b) induced by the definitions (* is the 3-Hodge star induced 
by the metric) 

F : = (*e) A dx° + b 
A : = aoda; + a. 

From here we build the hamiltonian version of the equations of motion; by using 
the rules 

* (* (dx 1 A dx j ) A da; ) = —dx i A dx 3 

* (dx 1 A dx 3 ) = * (da; 1 A dx 3 ) A da; 
the constraint submanifold will be given by the 3-integral sections of the EDS 

Jx° '■ = i* x ° (J\ ^x°) 

= (da-6,de = 0,d6 = 0) alg 

which are the usual constraints in the Hamiltonian version of electromagnetic field 
equations. 

5.3. Gotay-Nester algorithm. Before starting with the next examples, it is nec- 
essary to summarize the Gotay-Nester algorithm. The initial data is a triple 
(Mo, wo, fZo) where {Mq,ujq) is a presymplectic manifold and Hq G C°° (Mo). Then 
the algorithm proceed with the following steps: 

1. Calculate Kq := kerwo- 

2. For each I e N, let Ci +1 be the ideal in C°° (Mi) generated by Ci+i := 
{X-Hr.Xe Ki}. 

3. Define M; + i as the zero locus for functions in C;+i. Admitting that M/+i 
is a submanifold, define lu 1+ i := ui\ Ml+1 , H t+1 := Hi\ Ml+1 , and K t+1 := 

4. The algorithm stops whenever M; + i = Mi. 

5 That this system is linear is not a general fact; instead it is a property of the EDS which we 
are handled, and is equivalent to the linearity of the Maxwell equations. 
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5. The dynamics on the final constraint manifold Mi is governed by the solutions 
of Xh jwq = di^o tangent to Mi . 

Some shortcuts are in order here: Because M; + i C Mi, we have that (TMi) C 
(TMi + i) ± , so in each step there are dim (TMi + i)' L — dim (TMf) complementary 
vectors to look for. In order to find them out, it is useful the following remark: If 
F G Cj+i, and it admits a hamiltonian vector, then Xp G (TMj+i)^. 

5.4. Poisson sigma models. Here we will work in the setting of example 10. This 
example has the following feature, which is not shared by the rest of the examples 
we will deal with: Its Hamilton-Cartan EDS has torsion. In fact, it will be shown 
that the torsion is equal to the jacobiator of the Poisson structure tt. 

5.4.1. EDS analysis of Poisson sigma models. According to our previous discusion, 
we have that the Hamilton-Cartan EDS 1 is locally generated by 

(V := dx» - n^rjaud^, 

\T a := d Vaa A df* - | {d^ v ) v^rif^ A d£°. 

It can be shown (see for example [26]) that if J (tt) is the jacobiator of the bivector 
7r, then 

dr = [J (Tr)r p i la „ mp dC A d^ mod X. 

So the jacobiator in these systems is related to the torsion of the associated EDS. 
Because we will study the Cartan characters of a 2-dimensional integral element, it 
is not necessary to find an expression for dr M . 

Let us consider the flag E = C E\ := C E := (vi,v ), where 

Vi := d v + X?d xlt + ^ av d Vall G i = 0, 1. 

Then E will belong to V 2 (X) iff 

'X p a - n^vou = 0, 
XI - TT^Vlu = 0, 
[S? ff -Sj ff - (9^)770^ = 

for p, a = 1, • ■ • , dimAf ; then codim^V^ (I) = 3n. On the other side, for the polar 
spaces we have that 

H {E ) = {v G T (tx . v) S : vj (dx>* - ir^r/^dC) = 0, fx = I, ■ ■ ■ , dimAf} 
H (Si) = LeH (E ) : iuvu (d Va(7 A d£ a - ± (d a ^ v ) n^i^df A d£ Q J = 



|ij G H {E ) : di] la + terms in d^j, 



so the reduced characters for E will be Co (E) = dimM, c\ (E) = 2n. Then the 
Cartan test is satisfied, and the system is involutive at E. 

5.4.2. Investigation through the Gotay-Nester algorithm. It is interesting to see how 
the Jacobi identity for tt can be obtained from the Gotay-Nester algorithm associ- 
ated to the non standard variational problem describing a Poisson sigma model. So 
let us fix on £ a 1-submanifold Eq, which is locally described as £° = 0, and define 
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So := S\ So- According to our recipe, on T (So) there exists a presymplectic form 
and a hamiltonian; in order to get them, let us calculate 

dA = dr la ^dC Adx^ -(d^)r, a , lVp ,dx a AdC Ad^ -2n^ VPv d Va ^AdC Adf 



d\ 



So 



= -d?/^ A dx^ A d^ 1 



and 



<%>jA = rjo^dx 11 - Tr^rjo^rj^d^ d^j~X = i^dx^ - n^Tjo^d^ 1 . 

Sa 

In order to simplify notation, let us introduce the definitions r := £° , a := £ , r)n :— 
?]o^ and 7^ := r/i^; then if X, := (Sxf , 5rf v , Sj l p ) ,i = 1,2 indicates arbitrary vectors 
in T(,j. ,j 7 )T (Sq) the presymplectic structure reads 



(Xi,X 2 ) 
and the hamiltonian will be 



(5x^1 - 8x^1) da 



Ho ~ [ ^ [(xn' ~ ^ v lu] da. 



Here we denote with a prime the derivative respect to the cr-variable. So we have 
that 

kerw = {(0,^,0)} , 
and the primary constraints will arise from 

= dH (SVu) 

= f % [(O' - ir^jv] da => (x»)' - 7r^7, = 0. 

-'So 

Termination of the algorithm. Let C\ be the set C\ := {(a;'') - 71-^7,,} and M x C 
r (So) its zero locus. Now let us define 



where /„, \x = 1, • • • ,n are arbitrary functions on Eq and n = dim M; then 



dF (X 2 



S„ 



(SxZ)' -(d^^jJxZ-n^S^ 



da 



s„ 



and this function will have hamiltonian vector field, namely 
X F := (-/„*"", 0, (/ P )' + /ju (0,7^) 7 .) 
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The different choices for the arbitrary functions / M will gives the vector fields com- 
plementary to ker wo ■ The secondary constraints are obtained from 

= dH (X F ) 

= f V» [~ W)' + (dp^n l»f^ p - tt^ ((/„)' + U (d^) lp )] da 

= [ ^ [u (d P ir v n {xpy + (d p wn lu f u ^ - lp \ da 

J So 

- / »fo [f» ( d P^) + ( V) lvf^ P ~ ^"A. {duTf") 7p ] da 

JS 

= [ ri^fou, [(dp^n if + (d p n^) - tt"" (dp^)] da 

J So 

where was used that we are restricting on Mi, so that (x p )' = n^j^. Because 7r is a 
Poisson structure, these constraints are identically satisfied, and M\ is an invariant 
manifold. 

5.5. Dirac method and integrability conditions. From [29] we know that the 
following PDE system 

| ^zz + v4>xx = o 

\<j>yy = 

has the equations 4> xxy = <f> X xxx = as integrability conditions. This system can be 
written as the EDS I generated by 

|# := dip — pdx — qdy — rdz, 

T x := dr A dx A dy + ydp A dy A dz, T 2 := dg A dx A dzj 

on R with coordinates (x,y,z,(f),p,q,r). It is noted also that these integrability 
conditions can be obtained applying the prolongation procedure twice. In the fol- 
lowing sections we will use the correspondence settled by Theorem 4.5 in order to 
calculate the aforementioned integrability conditions through the Gotay-Nester al- 
gorithm. The idea is that, although the system initially considered is highly non 
regular, there exists (via Cartan-Kuranishi) a prolongation which is regular 6 , and 
the Dirac constraints can be determined by our theorem 4.5: Therefore if we know 
the Dirac constraints, we know something about the prolonged EDS. In particular, 
the integrability conditions must arise as these kind of constraints. 
Therefore we need to formulate the EDS as a non standard variational problem; 
thus let us define the double fibration 

A := R 7 Ai := R 4 M := R 3 

(x,y,z,(/),p,q,r) (->■ (x,y,z,4>) (->■ (x,y,z), 

taking as prolongations structure the given EDS I, and as lagrangian the trivial one 
A = 0. In order to build an associated Lepagean equivalent problem, we consider 

2 

A := A f\ (T*M) © R 2 



° There remains the question of regularity of the foliation introduced in order to define the 
presymplectic space; it is supposed that this condition is fulfilled. 
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with coordinates (x, y, z, <f>,p, q, r; a, A, /tt), a = AdxAdy+BdxAdz+CdyAdz, /i,A£ 
R; the functional to be extremized will be 



S := 



[ (a A 9 + ATi + fJ,T 2 ) • 
Jm 



' M 

We know that it is a covariant Lepage equivalent problem; the following lemma 
shows that it is also contravariant. 

Lemma 5.1. The projection 

v : A -> A : (x,y,z,(j),p,q,r;a,X,fi) ^ (x,y, z,<f>,p,q,r) 

maps extremals of S into integral sections for I; conversely, for each integral section 
of I given by 

(x,y,z) h-> (x,y,z,(j),p,q,r) 

there exists a £ ft 2 (M) , A, fx € C°° (M) such that joined to the integral section 
gives rise to an extremal of S. 

Proof. The variations respect to a, /i and A gives us the generators of X, so its 
projection will be an integral section of X. The variations respect to </>, p, q and r 
leads to 

'da = 0, 

ydA A dy A dz + a A dx = 0, 
d/i A d.T A dz + a A dy = 0, 
_ dA A dx A dy + a A dz = 0. 
This system is equivalent to 




in order to find a solution to it, we take A arbitrary, and define C and A by using the 
second and fourth equation respectively. From third equation we can determine B 
once fi is known, and the first of them gives us a differential equation for /i, namely 

I'-yy = + 2/A xx . 

So the system has solutions, and our new variational problem is contravariant. □ 

5.5.1. Associated dynamical system. Let A be the form below the integral sign in 
S, and let us define Nq C A as the submanifold consisting of the points in A such 
that z = 0. Then 

dA = da A (dcj) — pdx — qdy — rdz) — a A (dp A dx + dq A dy + dr A dz) + 

+ dX A (dr A dx A dy + ydp A dy A dz) + d^i A dq A dx A dz 

and thus 



= dA A dx A dy A dip + dA A dr A dx A dy; 

N 



dA 

this allows us to define on A^o the presymplectic form 

w z (Xi, X 2 ) := I [{6A x 8fa - 5A 2 5(j) 1 ) + {SXi6r 2 - 6X 2 6n)] dx A dy, 
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where we put 

Xi = (6(f>i, 8pi, Sqt, Sn; 5ai, 8Xi, SfXi) , i = 1,2 

in order to denote generic vectors in the tangent space of T (No). It is important to 
know the kernel of this form. 

Lemma 5.2. The kernel of ' uj z is given by 

keroj z = ((0, 8p, Sq, 0; 6Bdx A dz + SCdy A dz, 0, Sp,)) . 

Let us now build the hamiltonian; contracting A in the ^-direction 

d z j A = — (Bdx + Cdz) A (d<f) — pdx — qdy — rdz) — ra + y\dp A dy + fidq A da;, 

so that 

<9 z jA = 

No 

= - (Bdx + Cdy) A d(f) + (Bq - Cp) dxAdy- rAdx Ady + yXdp Ady + pdq A dx 
we obtain the following expression for the hamiltonian 



H z = / [(C4> x - B4> y ) + (Bq ~Cp- rA) + (y\ Px - fiq y )} dx A dy. 
Jr 2 

Primary constraints. The first order constraints can be obtained from 

dH z (keiuj 2 ) = 0. 
Taking into account the previous calculations, it results that 
= dH z (Sp) 

= / [yX (5p) x - C8p] dx A dy 
Jwl 2 

= (yK + C) 5pdx Ady => yX x + C = 

Jm. 2 

= dH z (Sq) 

dx A dy 

=> B + fly=0 



BSq - p(5q) y 

= (B + [i y ) Sqdx A dy 
Jr 2 

= dH z (SB) 

= / (q - 4> y )SBdx Ady =S> dj> y - q = 

Jr 2 

= dH z (SC) 

= (<t>x~ P) SCdx Ady => (f> x - p = 

Jr 2 

= dH z (Sp) 

q y Spdx Ady =>■ q y = 0. 



The set 

Ci := {y\ x + C,B + n y , <f) y -q,4> x - p, q y } 
gives rise to the primary constraints submanifold Mi = {C\ = 0} C T (No). In order 
to find out the elements in (TMi) 1 , let us use the remark made in the subsection 
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5.3: If F is a constraint, and it admits a hamiltonian vector, then Xp G (TMi)~ L . 
Let us take the constraints q y = y <p y — q = 0; none of them has a hamiltonian 
vector field, although the consequence (f> yy = admits such a vector. Therefore the 
function 

F\ -= fiwfidx Ady 



yyj 

JR 2 

(where f\ € C°° (M 2 ) is arbitrary) has a hamiltonian vector field. In fact, from the 
equation 

uj z (X Fl ,X 2 ) = dF 1 (X 2 ) 

and taking into account that 

dF 1 (X 2 )= [ (S<j> 2 ) hdxAdy 

JR 2 yy 

= [ S<f> 2 (h) dxAdy 

we conclude that 

X Fl = (o,0 ) 0,0;(/i) vy dsAdy,0 J o). 

Therefore we see that (TMi) ± = kerw z (B (Xp^, and we are ready to calculate the 
secondary constraints. 

Secondary constraints. These constraints are obtained from dH z uTJWi) 1 ! = 
and we know that on Mi, dfP (keru; z ) is identically zero; then 

= dH z (X Fl ) 



/ rih) dxAdy 



<- yy fidx Ady => r yy = 0. 
Let us now define C 2 ■= C\ U{r yy } and M 2 := {C 2 = 0}. By using the additional 

F 2 := / r yy f 2 dx A dy 
dF 2 (X 2 )= f (Sr 2 ) yy f 2 dxAdy= f 5r 2 (f 2 ) yy dx A dy; 

J R 2 J R 2 



constraint r yy we can define 



and so 



the associated hamiltonian vector field is 

X F2 = (0,0,0,0;0,(/ 2 ) w ,0). 

It is the complementary vector we were looking for, and so (TM 2 ) ± = (TMi)" 
(Xp 2 ). Thus we can calculate the new constraint, namely 

= dH z (Xp 2 ) 

= / (f 2 ) yp x dx Ady 
Jr 2 

= / h (yPx) vv dx A dy => (yp x ) 
Jr 2 
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But on M%, p = 4> Xl and therefore 

{yPx)yy = (y<f>Xx)yy 

= (y4>xxy + 4>xx) y 

— U^xxyy © ^tfixxy 

— ^(frxxy 

on M2. 

Tertiary constraints. Now C3 := C2 U {<f>xxy} and let M3 be its zero locus. Contin- 
uing with the algorithm, we need to find the vectors in the symplectic orthogonal 
to TM3 associated to the (new) constraints functions. So we define 

^3 := / f3<Pxxyd% A dy, 
Jr 2 

and its derivative along X2 will read 

diM*2)= I /;;: v d.r ' dy 

JR 2 

= -/ (h)xxy dfadx A dy. 

Then 

X F3 = (0, 0, 0, 0; - (f 3 ) xxy dx A dy, 0, 0) , 

therefore (TM3) = (TM2) © (Xf 3 ), and from here we obtain the new constraint 
= dH z (X Fa ) 

= / r(h) dxhdy 

JR 2 

= r xxy f 3 dx Ady => r xxy = 0. 

Jr 2 

Quaternary constraints. We define C4 := C3 U {r xxy } and M4 := {C4 = 0}. Then if 

Fa := / fir XX ydx A dy 
Jr 2 

we have that di<4 (X2) = — / R2 (fi) xxy Sr2dx A dy, and the new vector in (TM4) 
will be 

X Fa = (0,0,0, 0;0,- (h) xxy ,o) . 
Then we have that (TM4) 1 = (TM3) © {Xp 3 ), and for the new constraint 
= dff 2 (X F J 

= / (/i) IXS A dy 
Jr 2 y 

= (yp*)xxyU dx Ad v =>■ ,: .'//'•• : ',v ;/ °- 

But on M 4 

iyVx) xxy — (y& xx )xxy — (jj&xxxx) y — 4*xxxx © y^xxy — 4*xxxx: 

giving rise to the constraint 4> xxxx = 0. 
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Fifth order constraints. With the definitions C5 := C4 U {(f> xxxx } , M5 := {C5 = 0} 
and 

F5 ■= / f$4>xxxxd.x A dy 
Jm. 2 

we have that dH z (X 2 ) = J R2 {h) xxxx b&dx A dy, so 

X F5 = (0, 0, 0, 0; {h) xxxx dx A dy, 0, 0) . 

The new symplectic complement is given by (TM5) = (TM4,) 1 ' © (Xp 5 ). Thus it 
is obtained a new constraint from 

= dH»(X Fs ) 

t xxxx /s dx A dy =^ t xxxx 0. 

Termination of the algorithm. Let us define Mq C M5 as the zero locus of r xxxx 
in M5; by adding the hamiltonian vector field of Fq := J R2 fer xxxx dx A dy to 
(TMs) 1 , we obtain the vector space (TM§) J ~; it can be proved that the function 
dH z (Xp 6 ) is identically zero on M@, and we can conclude that it is an invariant 
submanifold of T (No). By performing an elimination of the auxiliar variables on 
the set of constraints Cq := C5 U {r xxxx } and keeping those involving only the 
original variables x, y, z, 4>, we see that the integrability constraints must be 4> xxy = 
0, 4> xxxx = 0. A note of caution about the applicability of the method is that this 
elimination process can in principle be as involved as the problem of finding the 
original integrability conditions. 

5.6. Non regular EDS and Gotay-Nester algorithm failure. We want to 
construct an example where the Dirac method fails in some sense; in order to do 
that, we will use the approach to variational problems developed here. The rationale 
behind this is the following: We know from the theorem 4.5 that whenever the EDS 
associated to the equations of motion admits regular elements, there exists a one- 
to-one correspondence between the solutions for this EDS and the solutions for the 
dynamical system on the sections of the selected spatial slice. Therefore we need 
to find a variational problem whose extremals are described for an EDS with non 
regular integral elements. Concretely, we will take a system where the non regularity 
stems from singularities of the set of integral elements: this lack of regularity cannot 
be circumvected via prolongation, and the Dirac method would fail. So the Gotay 
algorithm could run into problems in facing them. 

5.6.1. Initial setting. Let us consider the EDS I generated by 

6\ := du — qdx, 62 := dv — qdr 

on A := K 6 with coordinates (x, y, u, v, q, r); the independence condition is dxAdy ^ 
0. The idea is the one given in example 4, that is, by considering that X as the ideal 
defining a prolongation structure on A, which has the double fibration structure via 

A -> Ai := R 4 -> M := R 2 
(x,y,u,v,q,r) i-> (x,y,u,v) i-> (x,y), 

and by taking the lagrangian A = 0. Using our prescriptions, the Lepage equivalent 
will be constructed on A := (R 2 x M 4 ) © (T*M 2 © T*M 2 ), with action 

S:= f (aA0i+/3A0 2 ), 

JR 2 
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where a = adx + bdy, (3 = mdx + ndy. We need to prove that it is a bivariant 
Lepage equivalent problem. 

Lemma 5.3. The projection 

v : (x, y, u, v, q, r, a, (3) i-> (x, y, u, v, q, r) 

makes (K, 0, a A Q\ + f3 A 02^j into a Lepagean equivalent of (A, I, 0). Furthermore, 
it is bivariant, that is, every extremal of the variational problem posed by S on 
sections of A projects on an integral section of I (because the lagrangian is trivial, 
there are no action to work with here) and conversely, every integral section of I 
can be lifted to an extremal of S . 

Proof. The variations along a and (3 will give us the generators for I; this shows 
that it is a covariant Lepage equivalent problem. On the other side, variations along 
u,v,q and r leads to the following system on a and (3: 

da = 0, 
d/3 = 0, 

< 

(3 A dr + a A da; = 0, 
d(qp) = 0. 

For every (x,y) i— > (u,v,q,r) integral section for I, this systems admits solutions, 
so it is contravariant also. □ 



5.6.2. Gotay-Nester algorithm. Now we start with the analysis of the dynamical 
equations associated to the regular slicing y = 0. In order to do that, we take 
iVo C A the submanifold composed by the points in A such that y = 0; it is a bundle 
on M, and our next task is to define on its space of sections a presymplectic structure 
and a hamiltonian function. The first thing to note is that for A := a A 9\ + f3 A #2 

dA = da A (du - qdx) + a A dq A dx + d(3 A (dv - qdr) + [3 A dq A dr, 



so 

dA 

Therefore if 



N 



da A dx A du + dm A da; A (d?; — qdr) + mdx A do A dr. 



Xi = (Sui, 5vi,dqi, Sr^Sai, oft) , i = 1,2 
denotes a pair of arbitrary elements in TV (No), the presymplectic structure will be 



uy (x u x 2 ) 



Sa\Su2 — 8m\bvi — m8ri6q2 + 

+ (qSmi + mSqi) Sr 2 + SuiSa2 + (Svi — qSri) Sm 2 



dx 



and 



H v = / [b (u x - q) + n (v x - qr x )} dx. 
Jr 
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Primary constraints. As above 

kerw y = /(0,0,0,0;tf&dy,0),(0,0,0,0;0,<5ndy)\ 

and therefore the primary constraints will arise from dH v (keroj v ) = 0. Then 
= dH y (Zi) 

— J Sb (u x ~ q)dx => 

= dfl* (Z 2 ) 



Ux — q = 0, 



= / 5n (v x - qr x ) dx 



v x ~ qr x = 0. 



Let us define C\ := {u x — q, v x — qr x } and Mi := {C\ = 0}; in order to find the 
two vectors complementary to kerw y , we define 



fi (u x - q) dx, 



F-, 



h {v x ~ qr x ) dx; 



then 



X 2 -F x = J fi[(Su 2 ) x -8q 2 ] dx 

= - f [(fl) x Su 2 + f 1 Sq 2 ]dx, 

X 2 -F 2 = J f 2 [(Sv 2 ) x - r x Sq 2 - q (6r 2 )J dx 

= / [(9/2)2, Sr 2 - (/2)aj Sv 2 ~ r x f 2 6q 2 ] dx. 
Then, under the assumption m/0, we obtain that 

X Fl = fo,^,0 > £;(/ 1 ) SB d*,0 
\ m m 



., qr x h Ixh r x f 2 . . 
0, , , ;0, (f 2 ) x dx 



m m m 
yielding to the secondary constraints 
= dH v (X Fl ) 



qfi 
m 

Qxfi 



q 



di- 



rt) 

Q = dH y (X F2 ) 



d.r 



q x „ 
n — = 
m 



- b 



qxh 



qr x f 2 



q x f 2 

r x q 



Txf 2 



)]}■ 

/ xJ 



, qxh , 





b— = 0. 

m 



We are now in a very uncomfortable situation: Let us suppose that we define 
C 2 := CiU{nff,6^} and M 2 := {C 2 = 0}, and we want to calculate keruj y nTM 2 . 
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This set is composed by those vectors in kerw y that annihilates the differentials of 
the new constraints; therefore Z = (0, 0, 0, 0; 5bdy, Sndy) £ TM2 iff 

5n^ =0 

m 

5b^ = 0. 

m 

Thus if q x ^ the algorithm must stops, because in this case TM% = TM X , and 
it will continue in those points where q x = 0. This is a very bad behaviour of 
the Gotay-Nester algorithm, and it has to do with the singular nature of the new 
constraints. 



5.6.3. EDS analysis. The analysis of this PDE system is interesting, because its set 
of integral elements consists of several parts, each of them with different regularity 
behaviour. In fact, we will see that these differences introduces the singularities in 
the Gotay-Nester algorithm. 

The integral elements for X divides into two sets: the set W consisting of the 
elements with q x = 0, and Us is the set of integral elements such that n = b = 0. 
None of them has regular elements; in particular, the elements in U3 has 0-forms 
in its definition, so this set needs to be redefined. The procedure is to pull back 
the original EDS I to the submanifold defined by the collection of problematic 0- 
forms. It gives rise to a new EDS I3 whose set of integral elements admits a new 
decomposition into a pair of subsets V\ and V3. The set V3 is characterized by the 
0-form constraint m = 0, and will not be studied; the set V\ is composed by regular 
elements. Additionally the prolongation of the portion of the EDS corresponding 
to the elements in W gives rise to a new EDS Iw on G 2 (TR 10 ) which is regular 
too. 

Therefore it is necessary to know the main characteristics of our EDS X; first of all, 
a set of algebraic generators will be 

{61,62, dq A dx, dq A dr, da,d/3,/3 A dr + a A dx, /3 A dq} . 

Now we proceed to calculate the (reduced) characters for the flag Eq := C Et := 
(vi) C E := (vi,V2), where 

vi := d x + Uidu + Vtd v + Qid 9 + R^. + A x d a + B x d h + Mxd m + 

v 2 := d v + U 2 d u + V 2 d v + Q 2 d q + R 2 d r + A 2 d a + B 2 d b + M 2 d m + N 2 d n 

are in T( x y u >v q r )M. 10 . Then these components must satisfy 



Vij (du — qdx) = 






5 = 0, 


Uij (du — qdr) = 




Vi - 


qRl = 0, 


v 2 j (du — qdx) = 


<^> 


u 2 = 


0, 


v 2 j (dv — qdr) = 


44> 


v 2 - 


qR2 = 0, 


(vi A v 2 )jdq A dx = 




Q 2 = 


= 0, 


(vi A v 2 )jdq A dr = 




Q 2 R 


1 — QiR 2 


(vi A v 2 )_ida = 




A 2 - 


Bi=0, 


(vi A u 2 )jd/3 = 




M 2 - 


-N x = 0, 


(vi A v 2 )a (j3 A dr + a A dx) = 




nRi 


— mR 2 — 


A v 2 )jf3 Adq = 




nQi 


- mQ 2 = 
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'u 1 = 


q, U 2 = 0, 


Vi = 


qRi,V 2 = qR 2 , 


Qi = 


o,q 2 = o, 


A 2 - 


Si = 0, 


M 2 - 


- iVi = 0, 


nRi 


— mi?2 — 6 = 0, 



Thus V 2 (I) is described by the equations 

(U!=q,U 2 = 0, 
Vi = qRi,V 2 = qR 2 , 

Q 2 = 0,QiR 2 = 0, 
A 2 -B x = 0, 
M 2 - Ni = 0, 
ni?! — mi? 2 — & = 0, 
.nQ! = 0. 

It is necessary to analyze carefully this system. It results that V 2 (I) = W U f/ 3 , 
where 

(U 1 = q,U 2 = 0, 
V 1 = qR u V 2 = 0, 
Q 2 = 0,i? 2 = 0, 
( A 2 - B x = 0, 
M 2 -^i = 0, 
6 = 0, 
= 0. 

Let us calculate the polar spaces for the integral element E; we have that 
H ( E o) = {v £ T {x ^^v,q,r,a,b, m ,n)^- W ■ vj (du - qdx) = 0, vj (dv - qdr) = 0} , 
H (Ei) = ju e H (E ) : vj (vijdq A dx) = 0, vj (vi jdq A dr) = 0, vj (vijda) = 0, 
vj (vud/3) = 0, vj (du (P A dr + a A dx)) = 0, vj («i_i/3 A dg) = j 

= e If (So) : i'J (Qidx - dq) = 0, vj (Qidr - #idg) = 0, 

vj (Axdx + Bxdy - da) = 0, vj (Mi da; + A^dy - dm) = 0, 

vj (mdr - - bdy) = 0, vj (mdg - Qi(3) = o|. 

We can consider the reduced Cartan characters for E in each of the sets W and E/3; 
in every case we have that Co (E) = 2 and ci (S) = 6 on W, by supposing that 7 
ra/0. Moreover we have that codim^!!^ = 9, so there are no regular elements in 
W. 

On U3 we must take care of the fact that it is defined through 0-forms, namely, 
the constraints n = and b = 0. The regularity of I must be studied through the 
regularity of the EDS I3 obtained from I by pulling back to the submanifold 

A 3 := {(x, y, u, v, q, r, a, 0, m, 0) : x, y, u, v, q, r, a, m e K} ; 

a set of generators for I3 will be 

{81, 9 2l dq A d.T, dq A dr, da A dx, dm A dx, mdx A dr, mdx A do} . 

Then if 

W! := d x + UA + V x d v + Q x d q + R x d r + A x d a + M x d m , 
w 2 := d y + U 2 d u + v 2 d v + Q 2 d q + R 2 d r + A 2 d a + M 2 d m 



7 In points where m = 0, some of the characters decrease by one. 
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the set of 2-integral elements V 2 (I3) is described by the equations 



Ui =q,U 2 = 0, 
V 1 =qR 1 ,V 2 =0, 
Q 2 =0,Q 1 R 2 =0, 
A 2 = 0, M 2 = 0, 
mR 2 = 0. 



Thus we will have that 17$ = Vi U V3 , where 



^1 = < 



C/i = 9, 17a = 0, 

Q 2 = 0, R 2 = 0, 
U2 = 0,M 2 = 0, 



U x =q,U 2 = 0, 
V 1 =qR 1 ,V 2 =0, 
Q2 = 0,Qi = 0, 
A 2 = 0, M 2 = 0, 
m = 



(21) 



and codim^Vi = 8. On the other side, the polar spaces for the flag E' := C E[ := 
(wi) C E' := (wi,w 2 ) will be 

H {E ) = T^y^v^a,™,)^ 8 ■ Vj (du - qdx) = 0, vj (dv - qdr) = 0} , 
H(Ei) = 6 H (E ) : vj (wudq A das) = 0, vj {wudq A dr) = 0, 
tij (wxjda A da;) = 0, vj (wudm A da;) = 0, 

vj (wij (mdx A dr)) = 0, vj (w\j (mdx A dq)) = j 

= ju e if (E ) : vj (Qid.T - dq) = 0, vj (Qidr - .Rida) = 0, 
wj (Aidx — da) = 0, vj (Mi da; — dm) = 0, 

vj (mdr — mRidx) = 0, uj (mdq — mQidx) = o|, 

and for E £ V\ we have that Co (E) = 2, ci (£7) = 6. Therefore, the integral elements 
in V\ are regular. 

Finally it is necessary to consider the prolongation of the set W, in order to see if 
some of its elements can be considered as the tangent space of some solution. The 
prolongation of an EDS 1 C 0* (M) is the pullback of the contact EDS on G 2 (TM) 
to the set V 2 (I). Let us introduce the coordinates 

(x, y, u, v, q, r, a, 6, to, n, u x , v x , q x , r x , a x , b x ,m x , n x ,u y , v y , q y , r y , a y , by, TTLy, fly) 



on G 2 



lcr 



; then W is the subset described by 



W : I 



Ux = 


q, u y = 




V x = 


qr x ,v y = 


qr y , 


q x = 


0,q v = 0, 




a y - 


b x =0, 




n x - 


m y = 0, 




nr x - 


- TTtVy — b 


= 0, 
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and the pullback Tw of the contact system to W will have the generators 





= du - 


qdx 


'dO u = 







9y 


= dv — 


q (r x dx + r y dy) 


d6 v = 


— q {dr x A da; 


+ dr y A dy) 




= dq 




dd q = 







9 r 

< 


= dr 


r x dx - r y dy 
a x dx - a y dy 


dO r = 


— dr x A da; — 


dr y A dy 


9a 


= da — 


dd a = 


— da x A da; — 


da y A dy 


Ob 


= db- 


a y dx — b y dy 


d9 b = 


— da y A da; — 


dby A dy 


6m 


:= dm 


— m x dx — m y dy 


de m = 


— dm x A dx 


- dm y A dy 


@n 


:= dn - 


- m y dx — n y dy 


d9 n = 


—drriy A da; - 


- dn y A dy. 



Let E be the flag Eq := C E\ := (vi) C E := (v 1 ,v 2 ), with basis 

V! := d x + UA + Vid v + Q 1 d q + + Aida + Bidb + Ahd m + NA+ 
+ Ufa. + Vfa. + Q\d q!c + Rfa, + Afa x + Bfa x + Mfa x + Nfa x + 
+ Ufa y + Vfa s + Q l y d qy + R l y d ry + Ald ay + B\d by + Mfa ly + Nfa y , 

v 2 := d v + U 2 d u + V 2 d v + Q 2 d q + R 2 d r + A 2 d a + B 2 8 b + M 2 8 m + N 2 d n + 
+ Ufa, + + Qfa x + Rfa x + Mil, + Bfa + Mfa x + Nfa x + 

+ Ufa + V y 2 d Vy + Qfa y + Rfa y + Afa + B 2 y d by + Mfa y + Nfa y ; 

then V 2 (Iw) is given by 

U\ = q,U 2 = 0, Vi = qr x ,V 2 = qr y ,Qi = 0,Q 2 = 0,Ri = r Xl R 2 = r y , 
Ai = a x ,A 2 = a y , Bi = a y ,B 2 = b y , Mi = m x , M 2 = m v ,N\ = m y ,N 2 = n y , 
R 2 X - R\ = 0,A 2 x -Al = 0, Bl - B\ = 0, M| - M y =0,N 2 - R l y = 0. 

This means that codhneV^ {lw) = 21. For E integral element the polar spaces are 
H (E ) = |i> e TW : vj9 u = 0, vj9 v = 0, vj9 q = 0, iu9 r = 0, 

vj9 a = 0, iu9 b = 0, vj9 m = 0, vj9 n = o| 
H (Ei) = {veH (E ) : vj (dr x +•■•)= 0, vj (da x +■■■)= 0, 

vj (da y + •••)= 0,Vj (dm x + ■ ■ ■ ) = 0, vj (dm y + • • • ) = o|, 

and the reduced Cartan characters results cq (E) = 8, c\ (E) = 13; therefore E is a 
regular integral element. 

5.6.4. Return to Gotay-Nester algorithm. Nevertheless, it can be instructive to 
desingularize these constraints by mimicking the finite dimensional procedure: That 
is. by considering that the sections fulfilling the problematic constraints are the 
union of the sections such that n = b = with the set of sections satisfying the 
single constraint q x = 0. Despite the validity of such a procedure, we consider the 
continuation of the algorithm in each of these sets. 

• Case n = b = 0. In this case C' 2 := C\ U {n, b] and if M 2 is the zero locus for 
C 2 . Therefore the algorithm must stops, namely, M' 2 is an invariant manifold. 
In order to find the dynamical equations, let us obtain the hamiltonian vector 
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field for H y on F(No), restricting it to the submanifold M 2 . We will have 
that 



Xo-HV = 



|<5fe 2 (u x - q) + Sn 2 (v x - qr x ) + 

+ b [{5u 2 ) x - Sq 2 ] + n [{5v 2 ) x - r x 5q 2 - q [5r 2 ) x ] jdx (22) 

and so X 2 ■ H y \ M , — 0; if the hamiltonian vector field for H v reads X^v = 
(u y , v y , q y ,T y , dyAx + b y dy , uiyAx + rLydy) , the tangent solution to A/^ of the 
Hamilton eqs must be 8 Xhv — and the dynamical equations reads 



ft 







« y - gr^ = 
9y = 

Cly = 

m y = 0. 



It can be compared with the regular EDS V\ described by Eq. 21: This 
branch of the regularization procedure finds the solutions associated to regular 
elements in V\. 

• Case q x = 0. Now C' 2 ' := d U {q x }, and M' 2 ' = {C' 2 ' = 0}. If we define 



G 



its derivative along X 2 results 



gq x dx 



X 2 -G= g (Sq 2 ) x dx 



g x 5q 2 dx. 



Therefore 



9x 



9x 



X G = (o, 

and the invariance condition will reads 
= X G ■ H y 



9^,0,^,0; 0,0 G (TM' 2 y 
m m 



(«-) -«(-) 

V mix \vn) c 



nq x 
g x dx 

m 



dx 

q.< 



— =o. 



This condition is fulfilled on M% , so it is an invariant manifold for the dynamics 
defined by uj v and H y . Thus let us calculate the hamiltonian vector field Xhv , 
and then restricts it to M 2 : by taking into account Eq. 22 we have that 



X 2 ■ H y 



M'/ 



b x 5u 2 - {b - nr x ) 5q 2 - n x Sv 2 {nq) x Sr 2 



dx. 



8 A11 the hamiltonian vector fields for H v has the form Xjjv = Mi-Zi + /J2^2, and the tangent 
condition requires (11=^2 = 0. 
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Therefore the equations of motion are 

' u, u = 



Vy ~ 


qr y 


q y = 





a y = 


b x 


by = 


8b 


m y = 


= n x 


n y = 


Sn. 



In this case the solutions of the dynamical equations are in a one-to-one corre- 
spondence to solutions associated to the regular elements in the prolongation 
of W; moreover, the additional step in the Gotay algorithm has to do with 
the extra prolongation performed. 

5.7. Some remarks. It must be noted that there exists a field theory whose Dirac 
constraints has similar behaviour to the constraints found in the last example: It is 
the system composed of a massless spin 1 charged particle and the electromagnetic 
field [14]. A detailed analysis of such a system from the point of view adopted in 
this article is being performed; the true nature of the problems encountered there 
is under investigation yet [7]. Nevertheless, it is interesting to note here that the 
non standard perspective permits us to analyse any EDS (in particular, any PDE 
system) using techniques of (pre)symplectic geometry. The example 5.5 is a sample 
in this direction. 

6. Conclusions. In this paper it was built an entire dynamical theory for non 
standard variational problems, at least for those admitting a Cartan form. We were 
able also to describe the constraint submanifold in every case where such a descrip- 
tion is possible. It was proved that this submanifold is the space of integral sections 
of an EDS which is associated to the EDS generated by the Hamilton-Cartan equa- 
tions and the slicing of the space-time. We saw that the Dirac constraints has two 
sources: They can arise from the differential closure of the Hamilton-Cartan EDS, 
although some others can appears because of the involutiveness requeriment, which 
is a consequence of the stability condition defining the constraint submanifold. Four 
relevant examples showing these aspects were discussed in some extent. 
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Appendix A. On (formal) variational calculus. The purpose of this appendix 
is to set some common facts about geometrical concepts involved in the description 
of spaces of sections and the variational calculus used throughout the work [15, 16]. 

A.l. Variational derivative. Given the bundle F — > B, let us consider on T (F) 
a monoparametric family at of sections such that do = er; for every b 6 B we can 
calculate the following vector 



V(b) := A [trt 
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which has the following properties: 

• It is a vertical vector, that is, it belongs to ker (t„). 

• The map b i— > V (b) covers a, that is, it makes the following diagram commu- 
tative 

VF 

v 

B F 

Therefore V is a section of the pullback bundle a* (VF). So we make the following 
definition. 

Definition A.l. The tangent space at a of the space of sections T (F) of a bundle 
F B is the space of sections T (a* (VF)). 

Note. We want to remark that the previous one is just a definition, and does not 
implies the existence of differential structure in particular cases. 

Now we proceed in the opposite direction: Let us suppose that we have an 
element V £ T (a* (VF)), and we need to build a monoparametric family of sections 
cr f £ T (F) such that 

• They satisfies the initial condition o~q = a. 

• The section V can be reconstructed through the formula 

V(b):=±[a t (b)} . 
at t=o 

This is the case when V admits an extension 9 V £ T (VF); in general it is a problem 
with no solution. Nevertheless there exists a particular case for which there always 
exists such an extension, namely whenever F is a vector bundle. 

Proposition 6. Let F — > B a vector bundle, a £ T (F) a section and V & 
r (a* (VF)) a tangent vector at a . Then there exists V £ T (VF) which extends V . 

Proof. Because F is a vector bundle, we have the identification V a n,) (F) ~ t _1 (b), 
and then we can consider V as a section for F. Then defining the family of sections 
o~t '. b i — y o~ (b) -(- tV (b) we obtain a vector field on F (the tangent vector field) 
extending the given vector. □ 

A functional on F is a function on T (F) (from a set-theoretical point of view). 
We say that a functional S on F is differentiable at a £ T (F) iff for every extensible 
V e T (a* (VF)) the number 

V\„-S:=±(S[<fH (a)}) 

exists and it is independent of the extension. Here <j>t : T (F) — > T (F) is the flux 
in r (F) defined through ipt '■ F — > F, ipo = id^ associated to some extension for V 
(that is, <frt (a) : b M> (ip t ° cr) (b)). With respect to this subclass of functionals we 
define the variational derivative. 



9 This means that there exist on F a vector field V whose restriction to the image of a is equal 
to V 
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Definition A. 2. The variational derivative of the differentiable functional S in 
a £ r (F) is the form SS [a] : T a T (F) — » R, given by 

SS [a] (V) = V\ a ■ S. 

Remark 2. The uniqueness of the derivative of a functional depends strongly on 
each case; it has to do with the choice of the extension of the vector V. 

A. 2. Variational derivative for special functionals. We want to study the 
variational derivative of some class of differentiable functionals on F, that is, those 
that can be given by the formula 

S x [(t] ■= I o* (A) 

J B 

for some A £ 0" (F) , n = dim B. 

Proposition 7. [13] The variational derivative for S\ at a is given by 

SS X [a] (V) = [ a* (CyX) , 

J B 

where V £ T (VF) is some extension for V. In particular, on these kind of func- 
tionals the variational derivative is unique. 

Appendix B. Some facts about exterior differential systems. In this section 
we introduce the basic concepts of the theory of EDS; the sources of the material 
presented here are [5, 22, 23, 25]. The problem is to find out conditions that ensures 
us the existence of submanifolds on which a certain set of forms vanish. Since any 
of such submanifold annihilates the algebraic ideal generated by these forms and its 
differentials, let us introduce the following definition. 

Definition B.l. An exterior differential system on the manifold M is an ideal 
X C f2* (M) closed with respect to the exterior differential operator. We say then 
that I is differentially closed. A manifold which annihilates all the elements in I is 
called integral submanifold of the EDS. 

It will be supposed that the EDS does not contain 0-forms; it is not a true 
restriction in cases where the zero locus of these 0-forms is a submanifold of M. In 
order to build up the integral submanifolds of a given EDS, we find out the planes 
(of the corresponding dimension) that annihilates I, asking about the conditions 
under which these planes are tangent to some submanifold. The following definition 
is motivated by these ideas. 

Definition B.2. An integral element of dimension k for an EDS I is a subspace 
E C T X M of dimension k such that a\E = for all a £ I. The set of fc-integral 
elements associated to a given T will be denoted by Vk {T), and it is naturally 
included in the manifold Gk (TM), the /c-grassmannian on M. 

The verification that a given subspace is integral for some EDS involves the 
resolution of a (mostly large) linear system, as the following fact shows. 

Lemma B.3. Let us define I k :=ln O fc (M); moreover, let us fix 
(a{,--- ,al ± ,al,--- ,al 2 ,--- ,a{,--- ,«£ p ) C 0* (M) 

a set of generators (in the algebraic sense) for the EDS I such that {a\, • • • , <x l k } C 
I 1 for all I. Then 
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• V k (X) \ x = {E G G fe {T X M) :a\E = Va£ I fe } . 

• The subspace E C T X M is r-integral for X if and only if = for all 
1 < I < k m and 1 < m < r. 

It is important to note that in order to verify integrability in terms of algebraic 
generators of an EDS, it is necessary to prove that the given subspace annihilates 
all the generators of degree less than or equal to its dimension. 

If the EDS X is defined on a bundle F — ->• B, we ask about integral submanifolds 
which are sections of the given bundle, called integral sections 10 ; if n = dimB, the 
procedure that allows to find them works by fixing a nonzero element Q S fl n (B) 
and looking for n-integral manifolds for X on which r* (f2) ^ 0. This discussion 
motivates the following definition. 

Definition B.4. An EDS with independence condition on a manifold M is a pair 
(X, f2) composed by an EDS X and a differential n-form Q. The integral elements 
(resp. submanifolds) of such an object are the integral elements (resp. submani- 
folds) of X such that fl ^ on them. 

The independence condition, although apparently innocuous, changes dramati- 
cally the analysis of an EDS, as we will see in the discussion on the Cartan-Kahler 
theorem. 

B.l. Kahler-ordinary integral elements. First we need to speak about the no- 
tions of regularity which are necessary in the hypothesis of the Cartan-Kahler the- 
orem. The first thing that we have to impose is the regularity condition for the 
integral elements in the differential geometric sense. 

Definition B.5. Let X C O* (M) be a EDS. An n-integral element E E V n (X) 
is Kahler-ordinary iff there exists a neighborhood E G U C G n (TM) such that 
U n V n (X) is a submanifold of G„ (TM). 

This is a typical requeriment which ensures (at least locally) a description of the 
set of ?i-integral elements as the level set of some differentiable functions. In general 
the sets V n (X) are algebraic subvarieties of G n (TM) and so there exists singular 
points. 

B.2. Polar spaces and Kahler-regularity. In order to build up integral sub- 
manifolds of dimension n from an integral submanifold of (n — I)-dimensional, we 
need to consider the possible directions in which it is possible. 

Definition B.6. Given an element E £ T4 (X), we will define its polar space as the 
subspace generated for all the integral elements of dimension k + 1 associated to X 
which contains E. The polar space for an integral element E will be denoted by 
H(E). 

The following lemma gives us a way in which we can calculate a polar space. 
Lemma B.7. If E C T X X and {vi, • • • , v^] is a basis of E, then 

H(E) = {weT x X:<t>( Vl ,--- ,v k ,w)=0, V^I H1 }. 



The set of integral sections for an EDS on a fiber bundle will be denoted by the symbol T (I). 
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From the previous lemma we get some intuition about how the polar spaces are 
determined. That is, supposing that X has generators like in lemma B.3 and that 
E = (vi, • • • , v„); the vector x £ T X X will belongs to H (E) iff 

'al (v ils x) = 0, 
a l (vji.Vfc.x) = 0, 

< 

a™ +1 (vi, • • ■ ,v n ,x) =0, 1 < a < k n+1 . 

So, instead of having a set of algebraic equations, we have a linear system determin- 
ing x, which depends on a flag in E; the next notion to be introduced assures the 
regularity of this system, in the sense that it deals with the invariance of its rank. 

Definition B.8. A Kahler-ordinary element E £ V n (I) is Kahler-regular if 

codimi? (E) = codimi? {E') 

for all E' belonging to a neighborhood of E in V n (I) C G n (TX). 

We need to develop some strategies in order to deal with this condition, because 
it is difficult to verify in practice. The main theoretical resource in this direction is 
the Cartan's Test for Involutivity. 

Theorem B.9 (Cartan's Test for involutivity). Let Ek, < k < n a flag of integral 
elements for I at x, and let us define 

cj. := codimH (Ek) , < k < n — 1. 

Then 

codim En V n (I) > c + • • • + c n _i, 

where the codimension of V n (I) at E n is measured in the algebraic sense: It is 
defined as the maximum number of smooth functions on G n (TX) which vanish on 
V n (I) and has independent differentials at E n . 

Moreover, V n (T) is smooth of codimension exactly Co + • • ■ + c„-i at E n if and 
only if the Ek are all Kahler-regular for < k < n — 1. 

B.3. Cartan-Kahler Theorem. Once we have introduced all the necessary in- 
gredients, we can formulate a version of the Cartan-Kahler Theorem [22, 23]. 

Theorem B.10 (Cartan-Kahler Theorem). Assume T is an analytic EDS on M 
and P n C M is an analytic submanifold which is Kahler-regular and such that, at 
each p £ P, H (T P P) has dimension n + r + l. Moreover, let us assume that R £ M 
is an analytic submanifold of codimension r such that P £ R and T P R intersects 
transversely with H (T p P) . Then for each p £ P there is a neighborhood U £ R of 
p and a unique analytic (n + 1) -dimensional integral manifold N £ U containing 
PDU. 

The main drawback of this theorem is that we need to work in the analytic 
setting; for the systems which we consider here this hypothesis is fulfilled, although 
it is important to remark that it is a strong condition to be required. 



l<^<fc 2 ,l<ii<n 

1 < m < fc 3, 1 < h < :h < n 



(23) 
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B.4. Cartan-Kuranishi Theorem. According to the Cartan's test for involu- 
tivity, we can apply the Cartan-Kahler theorem only in those cases in which the 
codimension of the set of integral elements has the appropiate value. It can happen 
that the integral manifolds we are looking for cannot be obtained by using this exis- 
tence theorem. For example, to ensure that the set of 3-integral sections for an EDS 
can be obtained by restriction of 4-integral sections to some particular submanifold, 
the hypothesis of the Cartan-Kahler theorem needs to be verified for each 3-integral 
section. These ideas motivate the following concept, central in the theory of PDEs. 

Definition B.ll. Let Icfl' (M) be an EDS on M, and E 6 V n (1). A flag 
C Ei C • • • C E n = E is a regular flag iff Ek is a regular fc-integral element of X 
for 1 < k < n — 1. 

The regular flags allows us to solve the PDE behind an EDS as a sucession of 
Cauchy-Kowalevky problems; we can introduce a fundamental concept related to 
the integrability of the underlying PDE. 

Definition B.12. An EDS is n-involutive iff every n-integral element is the termi- 
nus of a regular flag. 

The main consequence of this definition and the Cartan-Kahler theorem is the 
following proposition (see for example [">]). 

Proposition 8. // an EDS is n-involutive, through any n-integral element passes 
an integral manifold of dimension n. 

One may prove there exists an operation, called prolongation, which allows us to 
associate to every EDS I on a manifold M and every n G N a new EDS I' on a 
manifold M' such that 

• There exists a submersion p : M' — >• M, and 

• the n-integrals elements of I and I' are in one to one correspondence via the 
map p. 

The key point is that, under certain technical assumptions [">], after a finite number 
of prolongations one can obtain an involutive EDS: This fact is known as Cartan- 
Kuranishi theorem. In this work we are assuming that these technical requeriments 
are fulfilled elsewhere. 
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